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ABSTRACT 


The  work  presented  Volume  II  is  concerned  with  the  structural 
mechanics  and  synthesis  of  graphite -fiber  reinforced  composite  materials. 
Optimization  studies  (minimum  weight  design)  of  stiffened  cylinders  similar 
to  the  fuselage  component  were  carried  out  by  two  different  techniques.  The 
behavior  of  composite  structural  elements  was  investigated  in  plate -buckling 
and  post-buckling  studies,  buckling  of  stiffened  plates  with  cutouts,  and 
application  of  shell  theory  to  anisotropic  cylinders.  Micromechanics  studies 
of  fibrous  composites  included  development  of  a  photoelastic  technique  for 
analyzing  frozen  stresses,  a  discrete  element  microstress  analysis  of  uni¬ 
directional  fiber  composites,  an  application  of  the  theory  of  physically 
nonlinear  elastic  solids  to  composite  materials,  and  an  investigation  of  the 
multiple  circular  inclusion  problem  in  plane  elastostatics .  Failure  of  com¬ 
posite  structural  elements  was  investigated  through  studies  of  the  failure 
mechanism  for  off-axis  composites,  the  fracture  toughness  of  composites, 
the  effect  of  fatigue  and  sustained  loads  of  cross-plied  composites,  the  notch 
sensitivity  of  cross -plied  composites,  and  the  effect  of  a  variety  of  known 
intentionally  interjected  defects  on  the  fracture  strength.  Further  studies 
were  concerned  'with  multiaxial  stress  testing  of  composite  cylinders.  The 
effects  of  circular  and  square  cutouts  in  flat  panels  and  in  a  stiffened  panel, 
and  methods  of  reinforcement  of  cutouts  were  investigated  and  correlated 
with  predictions  based  on  discrete  element  analysis. 
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SECTION  I 


INTRODUCTION 


The  present  program  is  a  continuation  of  the  work  performed 
under  Contract  AF  33(615) -31 10  (see  References  1  through'6).  This  report 
covers  the  fifth  (and  final)  year's  woik  of  a  program  which  represented  a 
novel  approach  designed  to  fulfill  three  different,  but  clearly  interdependent, 
needs  of  the  Department  of  Defense:  a  materials  need,  a  structural  design 
capability  need,  and  a  need  for  more  scientists  and  engineers  trained  in 
applied  materials  problems  and  advanced  design  methods.  The  Carbon 
Products  Division  of  Union  Carbide  Corporation,  Case  Western  Reserve 
University,  and  Bell  Aerospace  Company  have  formed  an  Association  to  meet 
these  needs. 


The  Association  has  formulated  a  broad  program  which  includes 
the  development  of  new  materials,  generation  of  advanced  analyses  and  design 
methods,  and  education  of  graduate  students.  In  brief,  the  major  objectives 
are  (1)  to  develop  high  modulus  graphite  fiber  composites,  (2)  to  extend  the 
methods  of  structural  mechanics,  (3)  to  identify  DOD  applications  toward  which 
the  program  efforts  should  be  directed,  (4)  to  educate  engineers  capable  of 
developing  and  using  modern  materials,  and  (5)  to  integrate  materials  research 
with  the  needs  of  the  design  by  extending  the  technique  of  structural  synthesis 
to  include  material  variables. 

The  primary  responsibilities  of  Union  Carbide  Corporation, 

Carbon  Products  Division,  are  the  development  and  production  of  composite 
materials  and  the  measurement  of  those  mechanical  and  thermal  properties 
needed  for  the  structural  design  work  within  the  Association.  The  technical 
program  at  Union  Carbide  consists  of:  (1)  materials  research5  a  basic  research 
program  to  develop  new,  improved  composites  of  high  modulus  graphite  fibers 
in  both  resin  and  metal  matrices;  (2)  materials  fabrication,  an  applied  research 
program  to  produce  materials  for  the  joint  research  programs  of  the  Associa¬ 
tion  and  to  seek  new  ways  of  fabricating  components  which  better  utilize  the 
superior  properties  of  composite  materials;  (3)  properties  evaluation,  the 
measurement  of  the  mechanical  and  thermal  properties  of  certain  composites 
to  provide  data  for  the  joint  research  programs  of  the  Association;  and  (4)  fail¬ 
ure  criteria,  a  basic  research  program  to  determine  experimentally  adequate 
failure  criteria  for  anisotropic  materials  under  multiaxial  stress  states  and 
to  find  ways  of  representing  the  failure  surface  which  can  be  used  by  the  designer 
in  practical  calculation. 

The  work  at  Case  Western  Reserve  University  has  two  major 
objective i.  The  first  objective  is  to  advance  the  basic  structural  mechanics 
technology  required  for  rational  design  with  composite  materials.  Composite 
materials  offer  the  structural  design  engineer  the  prospect  of  being  able  ulti¬ 
mately  to  carry  on  simultaneously  the  design  of  the  structural  configuration 
and  the  material.  Achieving  this  capability  will  require  fundamental  advances 
in  structural  synthesis  as  well  as  a  substantially  improved  understanding  of 
the  behavior  of  composite  materials.  The  goals  of  the  structural  mechanics 
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research  program  at  Case  are  (1)  the  quantitative  formulation  and  efficient 
solution  of  the  structural  synthesis  problem,  including  material  variables, 
for  elementary,  but  representative  components  fabricated  from  composite 
materials;  (2)  experimental  stress  analysis  studies  and  theoretical  investi¬ 
gations  in  micromechanics  with  the  objective  of  improving  the  measurement 
and  calculation  of  stiffness  properties  and  failure  mode  criteria  for  compos¬ 
ite  materials;  and  (3)  the  development  of  improved  analysis  methods  for 
anisotropic,  nonlinear,  and  nonconservative  materials.  The  second  objective 
of  the  work  at  Case  Western  Reserve  University  is  to  develop  new  or  improved 
graphite  fiber-resin  composites  through  materials  research.  At  present,  the 
knowledge  of  fiber  surface  morphology  and  the  relation  between  fiber  surface 
characteristics  and  interfacial  adhesion  to  the  resin  system(s)  is  incomplete. 

A  better  understanding  of  these  interfacial  interactions  will  lead  to  improve¬ 
ments  in  presently  used  fiber-resin  composites  and  will  ultimately  permit  the 
judicious  selection  of  new  resins  and  new  fabrication  methods,  thus  leading  to 
a  second  generation  of  advanced  composites. 

The  primary  purposes  of  Bell  Aerospace  Company's  participa¬ 
tion  in  this  program  are  to  interject  user  requirements  into  the  applied  mate¬ 
rials  research  efforts;  to  apply  at  the  prototype  design  level,  the  advanced 
analytical  procedures  and  improved  understanding  of  material  behavior  which 
will  result  from  the  research;  at  e  to  establish  application-related  property 
specifications  for  materials  research  activities.  To  attain  these  objectives, 
a  six-part  technical  program  is  being  performed  by  Bell;  (1)  application 
selection,  the  objective  of  which  is  to  define  representative  configurations 
and  environmental  conditions  which  reflect  DOD  requirements;  (2)  recogni¬ 
tion  of  failure  modes,  a  task  which  involves  the  overall  structural  behavior 
such  as  elastic  instability,  deformation  limits,  and  fracture  and  the  material 
failure  modes;  (3)  determination  of  the  nature  of  and  methods  for  the  applica¬ 
tion  of  analytical  tools  needed  to  cope  with  the  anisotropic,  anelostic,  and 
nonconservative  material  property  behavior  and  the  multiaxial  stress  dis¬ 
tributions  anticipated  in  structural  configuration  associated  with  the  use  of 
the  subject  materials;  (4)  structural  synthesis,  a  task  which  involves  the 
application  of  structural  synthesis  techniques  at  the  practical  level  to  define 
the  most  desirable  material  compositions  within  a  particular  class  of  com¬ 
posites,  (5)  study  of  creative  design  concepts  which  will  be  required  because 
of  the  complex  material  behavior  of  composites;  (6)  testing  to  verify  the  value 
of  analysis  procedures  used  to  design  composite  materials  and  the  components 
made  from  composite  materials. 

Although  this  report  of  the  Association's  fifth  year  activities  is 
the  last  report  of  the  Association,  the  Case  Western  Reserve  University  will 
be  continuing  research  on  composite  materials  during  the  sixth  year  under 
sponsorship  of  the  Advanced  Research  Projects  Agencv. 

This  report  is  divided  into  three  volumes.  Volume  I  covers  the 
effort  on  materials  resea  rch;  Volume  II  coveis  the  work  on  structural  mech¬ 
anics,  analysis,  and  optimization;  and  Volume  III  covers  the  testing,  stress 
analysis, and  evaluation  of  a  representative  subscale  fuselage  component.  This 
division  was  made  because  particular  projects  might  be  of  interest  to  a  parti¬ 
cular  audience;  this  arrangement  also  made  possible  the  reduction  of  the  physi¬ 
cal  size  of  each  volume.  However,  the  following  Summary  (Section  II)  also 
contains  a  brief  outline  of  the  contents  of  the  other  two  volumes. 
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SECTION  II 


SUMMARY 


The  technique  of  structural  synthesis,  or  automated  optimum 
design  of  structural  systems,  has  been  extended  to  include  both  geometric 
and  material  parameters  in  the  minimum  weight  design  of  stiffened  fiber 
composite  cylindrical  shells.  A  special  computation  algorithm  was  devel¬ 
oped  for  dealing  with  design  parameters  such  as  ply  angles  which  do  not 
explicitly  appear  in  the  weight  function. 

Use  of  finite  element  methods  of  analysis  in  the  optimization  of 
composite  structures  was  found  a  costly  procedure.  Because  each  solution 
is  lengthy  and  numerous  solutions  are  required  during  the  optimization, 
this  approach  becomes  prohibitive.  An  alternate  method  of  analysis  based 
on  engineering  beam  theory  should  afford  a  more  economical  solution. 

A  Ritz  method  of  solution  has  been  applied  to  the  post-buckling 
analysis  of  laminated  fiber  composite  plates.  Complete  post-buckling 
response  curves  were  obtained  for  symmetrically  and  unsymmetrically 
laminated  plates  with  either  clamped  or  simply  supported  boundary  condi¬ 
tions.  The  results  are  used  io  define  an  effective  width  for  the  post-buckling 
strength  analysis  of  high  performance  composite  materials.  An  experimental 
program  to  verify  the  post-buckling  capacity  of  composite  plates  was  initiated. 
An  experimental  study  of  the  effects  of  cutouts  on  the  buckling  loads  of  lami¬ 
nated  plates  was  initiated  but  a  considerable  amount  of  time  was  spent  on  the 
fabrication  of  quality  specimens.  The  study  of  the  applicability  of  various 
shell  theories  to  anisotropic  fiber  composite  cylinders  has  been  completed. 

The  exact  solution  for  the  eigenvalue  buckling  of  anisotropic  shells  subjected 
to  axial  compression,  lateral  pressure,  or  gross  torsion  has  been  obtained 
and  programmed  for  the  computer. 

The  development  of  a  suitable  photoelastic  technique  for  the 
analysis  of  filamentary  composites  using  frozen  stresses  and  transmitted 
light  has  been  completed.  The  frozen  stress  technique  was  also  studied  for 
adaptation  to  scattered  light  photoelasticity  but  the  low  fringe  orders  of 
epoxy  matrices  necessitated  the  use  of  the  live-load  approach. 

A  special  finite  element  for  the  microstress  analysis  of  fila¬ 
mentary  composites  has  been  developed,  verified  on  problems  with  known 
solutions  and  applied  to  a  variety  of  common  problems  in  composite  mate¬ 
rials  for  which  solutions  do  not  exist.  The  results  of  this  study  are  used 
to  obtain  stress  and  strain  concentration  for  use  in  conjunction  with  a  mech¬ 
anistic  failure  theory  developed  under  this  program.  The  general  solution 
for  the  elastic  curvilinear  inclusion  problem  has  been  found  in  the  case  of 
antiplane  or  longitudinal  shear  deformation,  and  has  been  utilized  for  the 
analysis  of  randomly  spaced  fibers  in  a  filamentary  composite.  Plane 
deformation  of  a  similar  model  is  currently  under  investigation.  Pre¬ 
liminary  results  on  the  validity  of  classical  laminate  plate  and  shell  theories 
when  applied  to  fibrous  composites  have  been  obtained. 
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A  failure  mechanism  for  off-axis  fiber  composites  has  beer  postu¬ 
lated  based  on  extensive  experimental  observation.  It  was  observed  that  the 
matrix  is  in  the  bulk  plastic  region  early  in  the  loading,  consequently  failure 
is  hypothesized  as  a  critical  creep  rate.  The  state  of  stress  in  the  matrix  is 
calculated  using  the  concentration  factors  and  the  Neuber  leading  function 
analysis. 


The  fracture  toughness  of  unidirectional  fiber  glass  composites  has 
been  studied.  Its  dependence  on  specimen  width,  crack  depth,  and  strain 
rate  has  been  determined  experimentally.  The  dependence  of  fracture  tough¬ 
ness  on  fiber  orientation  for  crossplied  specimens  is  also  being  studied. 

Crack-free,  orthogonal,  seven-ply  "Thornel"  50S,  ERLA4617 
composites  were  tested  under  tension-tension  fatigue  and  sustained  load  con¬ 
ditions.  The  density  of  cracks  resulting  from  these  tests  was  determined 
as  a  function  of  load,  time-at-load  and  number  of  fatigue  cycles.  These 
cracks  were  found  to  be  similar  to  those  frequently  developed  due  to  thermal 
stresses  on  cooling  from  the  curing  or  pressing  temperature  and  appeared 
to  have  little  influence  on  the  modulus  or  ultimate  tensile  strength  of  the 
composites.  Crack  density  determinations  as  a  function  of  stress  during 
the  first  load  cycle  showed  the  existence  of  a  threshold  for  crack  creation  at 
a  strain  of  approximately  50  percent  of  the  ultimate  tensile  strain,  fatigue 
cycling  to  105  cycles  at  a  load  of  85  percent  UTS  produced  no  measurable 
degradation  in  tensile  strength,  yet  yielded  crack  densities  higher  than  those 
found  in  samples  tested  to  failure  in  a  single  cycle.  In  sustained  load  tests, 
cracking  was  found  to  vary  nearly  linearly  with  stress  in  samples  loaded 
for  1000  hours.  Samples  loaded  for  shorter  times  showed  that  the  crack 
density  doubles  for  each  increase  by  a  factor  of  10  in  time  at  load.  No 
significant  creep  occurred  during  the  sustained  load  tests. 

Comparative  notch  sensitivity  studies  were  carried  out  on  orthog¬ 
onal,  seven-ply  "Thornel"  50S,  ERLA  4617  epoxy  and  P-1700  polysulfone 
composites.  Since  polysulfone  resins  are  tougher  and  have  much  higher 
elongations  to  failure  than  the  relatively  brittle  epoxies,  it  was  expected 
that  the  polysulfone  composites  should  be  less  notch  sensitive.  The  experi¬ 
ments  confirmed  this  expectation. 

Tensile,  flexural,  and  compression  tests  were  performed  on  a 
number  of  orthogonal  "Thornel"  40,  ERL  2256  composites  containing  inten¬ 
tional  defects  such  as  mylar  and  yarn  inclusions  and  butt  or  lap  joints. 
Reductions  in  tensile  strength  were  noted,  while  the  effects  on  flexural  and 
compressive  strength  were  generally  minor. 

Fracture  strengths  under  uniaxial,  torsional,  and  combined 
stresses  have  been  determined  for  "Thornel"  50S/ERLB  4617/MPDA  com¬ 
posite  cylinders  with  a  four-ply  (90°  ,  0° ,  0° ,  90°  )  lay-up  pattern.  The 
uniaxial  tensiP  find  compressive  strengths  of  the  same  composite  were 
measured  using  flat  coupon  specimens,  while  shear  strength  was  meas¬ 
ured  by  the  torsion  bar  method.  The  fracture  strengths  under  combined 
stresses  were  measured  by  subjecting  cylindrical  specimens  to  combined 
loadings  of  internal  pressure,  axial  force,  and  torque.  These  data  establish 
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t-he  approximate  shape  for  the  fracture  strength  envelope  of  the  composite 
encompassing  the  shear-tensile-tensile  stress  state. 

Effects  of  circular  and  square  cut-outs  in  flat  panels  and  in  a 
stiffened  panel  assembly  were  investigated.  An  unreinforced  hole  in  a 
multidirectional  laminate  test  specimen  was  shown  to  reduce  load  capability 
by  a  factor  of  two.  Integral  composite  reinforcements  around  holes  in  sev¬ 
eral  tensile  specimens  were  found  to  efficiently  carry  loads,  i.  e.  ,  develop 
the  strength  of  the  panel  without  a  hole.  The  latter  specimens  were  lighter 
in  weight  than  comparable  strength  specimens  with  "bonded  on"  external 
reinforcements  or  integral  internal  metal  reinforcements.  Strain  responses 
and  ultimate  loads  compared  reasonably  well  with  predictions  based  on  dis¬ 
crete  element  analysis.  A  technique  of  lightweight  reinforcement  for  a  cut¬ 
out  in  a  stiffened  compression  panel  was  successfully  demonstrated  in  test. 

Contents  of  Volume  I; 

The  work  presented  in  Volume  I  is  concerned  with  materials 
research  on  graphite-fiber  reinforced  composites.  Graphite*fiber  surfaces 
were  characterized  by  gas  phase  and  solution  adsorption  experiments  and  by 
Raman  spectroscopy.  The  latter  technique  allows  a  differentiation  between 
fibers  of  different  origins  and  heat-treatment  temperatures.  The  fabrica¬ 
tion  of  "Thornel"  fiber,  polyamide -imide  and  polysulfone  composites  was 
investigated.  Fabrication  of  polyamide-imide  composites  was  very  tedious 
and  evaluation  was  limited  to  determinations  of  torsion  shear  strength. 
Several  plates  of  polysulfone  matrix  composites  were  fabricated;  the  eval¬ 
uation  of  these  plates  is  presented  in  Volume  II  of  this  report.  "Thornel" 
fiber  composite?  were  further  prepared  by  in  situ  polymerization  of  nylon. 
The  epitaxial  crystallization  of  nylon  monomer  on  graphite  fiber  surfaces 
has  been  investigated.  Graphite-fiber,  nickel-matrix  composites  were 
further  characterized  at  room  and  at  elevated  temperatures.  Measure¬ 
ments  of  Young's  moduli,  tensile  strength,  stress  rupture,  thermal  expan¬ 
sion,  and  thermal  fatigue  are  presented.  Graphite  filament  degradation 
resulting  from  fabrication  (hot-pressing)  of  nickel  matrix  composites  was 
also  determined. 

Contents  of  Volume  III: 


The  work  presented  in  this  volume  is  concerned  with  the  per¬ 
formance  prediction,  testing,  and  post-test  evaluation  of  a  representative 
graphite-fiber,  resin-matrix  aircraft  fuselage  component.  Additional  mate¬ 
rial  properties  were  determined  and  structural  margins  of  safety  defined  by 
discrete  element  analysis.  After  seven  response  tests  under  various  load 
combinations,  the  component  was  tested  to  destruction  under  combined 
bending  and  shear  loads.  Failure  occurred  at  110  percent  of  the  target 
design  load  but  below  the  failure  load  predicted  from  tests  on  flat  panels. 

A  weight  saving  of  27  percent  over  an  aluminum  structure  of  equivalent 
strength  was  demonstrated.  The  component  was  also  three  times  stiffer 
than  an  aluminum  structure  of  the  same  weight.  Performance  projections 
indicate  that  the  same  component  built  with  presently  available  "Thornel"  50S 
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fibers  would  offer  a  weight  saving  of  49  percent.  Post-teat  evaluations 
included  tensile  and  compression  tests  on  curved  skin  panels*  on  stringers 
and  on  stringer-skin  combinations.  Optical  and  electron  microscopic  exam 
ination  of  the  fracture  surfaces  provided  further  insight  into  the  failure 
mechanism. 
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SECTION  III 


COMPOSITE  SYNTHESIS 


A.  Stiffened  Cylinder  Optimization 

(Professor  kicher  and  br.  T.  L.  Chao) 

With  the  development  of  new  composite  material s>  a  desiqn  enqineer  has  th? 
opportunity  of  ultimately  being  able  to  design  a  structural  configuration  and 
material  simultaneously.  However,  to  accomplish  this  task  it  will  require 
fundamental ,  advances  in  structural  synthesis  as  well  as  an  improved  understand¬ 
ing  of  the  structural  behavior  of  composite  structures.  This  work  is  aimed  at 
advancing  the  techniques  of  structural  synthesis  in  design  with  fibrous  com¬ 
posites.  By  extending  the  technique  to  composition  and  material  variables  it 
provides  a  tool  to  assist  the  engineer  in  the  intelligent  exploitation  of  the 
potential  offered  by  composite  materials. 

c  ■  Jhe  b*sAc  concePts  of  structural  synthesis  were  developed  and  reported  by 
Schmit  and  his  associates.  (7,8)  The  concepts  have  been  applied  successfully 
to  the  optimum  design  of  integrally  stiffened  cylindrical  shells  by  Kicher(9) 
ana  Schmit-Morrow-Kicher.  (10)  In  general,  structural  synthesis,  or  auto¬ 
mated  optimum  design  of  structural  systems,  can  be  cast  into  a  model  of  a 
mathematical  programming  problem  which  seeks  to  maximize  or  minimize  a  func- 
tion  of  several  variables.  The  variables  are  subjected  to  certain  constraints 
which  are  not  amenable  to  solution  by  the  classical  methods  of  calculus. 
Composite  structures  must  be  designed  with  an  integrated  design-analysis 
approach  supported  closely  by  materials  and  processing  criteria.  In  order  to 
carry  out  this  work  successfully,  a  designer  is  required  to  have  a  comprehen¬ 
sive  knowledge  of  design,  manufacturing,  and  analytical  methods.  This  compli¬ 
cates  greatly  the  design  processes.  It  is  hoped,  in  this  respect,  that  an 
extended  technique  of  structural  synthesis  will  assist  an  engineer  to  achieve 
an  optimum  structural,  design. 

The  main  features  of  the  computer  capability  developed  for  findinn  a  mini¬ 
mum  weight  design  of  stiffened  fiber  composite  cylinders  can  be  summarized  in 
the  following: 


1.  Multiple  load  conditions  -  each  load  condition  composed 
of  axial,  radial,  and  torsional  load  combinations. 

2.  Design  begins  with  the  elastic  properties  of  the 
constituent  materials,  fiber  and  matrix. 

3.  Design  variables  include  the  material  parameters  of 
fiber  orientations  and  fiber  content. 

4.  The  buckling  analysis  of  the  heterogeneous  anisotropic 
cylindrical  shell  is  based  upon  the  small  deflection 
theory  of  thin  shells.  The  cylinder  is  assumed  to 
buckle  into  a  torsional  wave-form. 
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5.  Eight  sets  of  boundary  conditions  are  provided 

6.  Pounds  on  all  design  variables  are  permitted. 

7.  Optimum  design  is  obtained  by  minimizing  the  objective 
function  with  the  design  variables  subject  to  side  and 
behavior  constraints. 

Figure  1  shows  a  sketch  of  the  shell  with  stiffeners  located  inside  the 
cylinder  wall.  The  design  variables  are  the  configuration  parameters: 
dimensions  and  spacings  of  the  stiffeners;  and  the  material  parameters:  fiber 
orientations  and  fiber  content.  The  thickness  of  the  composite  is  determined 
from  the  fiber  volume  content.  The  stiffened  shell  is  designed  against  the 
following  modes  of  failure. 


1. 

gross  buckling 

2. 

panel  buckling 

3. 

skin  buckling 

4. 

stiffener  bucklings 

5. 

combined-stress  failure  in  a 

ply 

6. 

delamination  failure  between 

two  adjacent  plies. 

Behavior  constraints  are  formulated  based  on  the  failure  stresses  and  strains. 
Side  constraints  are  the  lower  and  upper  bounds  on  the  design  variables,  and 
the  limits  on  the  geometric  admissibilities  among  the  dimensions  and  spacings 
of  the  stiffeners.  Together  they  form  the  constraint  set.  The  objective  func¬ 
tion  is  the  total  weight  of  the  cylinder. 

If  the  stiffeners  in  both  the  longitudinal  and  circumferential  directions 
are  assumed  to  be  closely  spaced,  then  the  elastic  properties  of  the  stiffeners 
may  be  considered  to  be  smearable  in  the  two  directions  over  the  skin.  Thus 
the  force-displacement  relations  are  obtained  by  performing  the  appropriate 
integrations  of  the  stresses  over  the  skin  and  the  stiffeners  with  the  mid¬ 
surface  of  the  skin  taken  as  the  reference  surface.  Timoshenko's  buckling 
equilibrium  equations  (11)  are  used  for  the  buckling  analysis  of  the  cylinder. 
The  skin  buckling  loads  are  approximated  with  Donnell's  equations.  (12)  The 
hat  section  stiffeners  are  analyzed  as  a  collection  of  plate  elements.  Methods 
to  determine  the  limiting  strength  of  composite  materials  are  based  on  the  work 
of  Chamis.  (13)  Material  failure  modes  considered  are  the  individual  ply 
failures  and  the  ply  delaminations  between  two  adjacent  plies. 

The  minimum  weight  design  problem  is  cast  into  a  model  of  a  nonlinear 
mathematical  programming  problem.  The  Fiacco-McCormick  penalty  function  tech¬ 
nique  is  used  for  so'ving  the  programming  problem.  Since  the  objective  or 
weight  function  is  independent  of  the  ply  orientations,  special  modifications 
have  been  made  in  the  synthesis  scheme  so  that  the  weight  is  minimized  with 
respect  to  all  design  variables  and  the  strength  of  the  structure  is  maxi¬ 
mized  with  respect  to  the  ply  angles. 
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Figure  1.  Integrally  Stiffened  Cylinder  with  Hat  Cross-Section  Stiffeners. 
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A  topical  report  on  this  subject  has  been  completed  by  T.  L.  Chao  (14)  and 
submitted  to  Case  Western  Reserve  University  in  partial  fulfillment  of  the  Ph.D. 
Degree.  That  document  is  also  available  as  an  AFML  report  TR  69-251.  The 
results  were  summarized  and  presented  at  the  AIAA/ASME  11th  Structures,  Struc¬ 
tural  Dynamics  and  Materials  Conference,  Denver,  Colorado,  April  22-24,  1970. 
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B.  Minimum  Weight  Design  of  Composite  Structures 

(Dr.  R.A.  Geliatly  and  D.  Dupree/  Bell  Aerospace) 

Development  of  an  automated  optimization  program  for 
stiffened  composite  panels  was  undertaken.  The  basis  for  the 
program  was  utilization  of  discrete  element  structural  analysis 
methods  in  combination  with  a  mathematical  programming  penalty 
function  formulation  This  approach  differs  from  that  of  'the 
previous  report  section  which  considered  a  "smeared  stringer" 
model  of  the  fuselage  structure.  It  was  recognized  from  the 
outset  that  the  introduction  of  composite  materials  would  z,?-?«e 
the  level  of  complexity  considerably  over  an  equivalent 
optimization  program  operating  on  conventional  metallic 
structures . 

The  overall  program  logic  flow  was  developed  and 
incorporated  a  streamlined  version  of  an  operational  Bell  General 
Purpose  Analysis  Program  combined  with  an  unconstrained 
optimization  procedure.  Because  of  the  complexity  of  the 
resulting  large-scale  program,  development  and  check-out 
procedures  were  conducted  in  two  stages.  In  the  first  stage, 
attention  was  confined  to  use  of  simple  isotropic  material 
elements  to  verify  operation  of  the  optimization  procedure.  In 
the  second  stage,  analysis  procedures  were  expanded  to 
accommodate  anisotropic  material. 

A  sample  optimization  problem  was  attempted  for  the  first 
stage  check-out,  and  significant  weight  reductions  were  achieved 
when  geometric  variables  and  member  sizes  were  used  as  design 
parameters.  The  problem  consisted  of  a  two-dimensional  isotropic 
membrane  panel  fixed  along  one  edge  and  loaded  in  its  own  plane. 
The  manner  of  support  and  loading  was  designed  to  represent  a 
tip- loaded  cantilever  beam-plate.  Figure  2  shows  the  finite 
element  idealization  of  the  initial  design  which  consisted  of  12 
rectangular  membrane  plate  elements.  The  structure  was  required 
to  be  symmetric  about  its  horizontal  center  line,  and  analysis 
was  confined  to  the  upper  half  of  the  beam  with  appropriate 
boundary  conditions.  Design  variables  were  the  y-coordinates  of 
Nodes  1-8  and  two  thicknesses:  the  three  upper  elements  of  the 
top  beam  half  were  required  to  have  equ*l  thickness  ti  and  the 
three  lower  elements  were  required  to  have  equal  thickness  t2  . 
Figure  3  shows  the  design  achieved  after  three  minimizations1. 
The  initial  and  final  values  of  the  design  variables  are  given  ^.n 
Table  I.  The  corresponding  initial  and  final  weights  are 
3 s 030  and  0,216  lb,  respectively.  The  initial  design  does  not 
represent  an  efficient  starting  point  such  as  a  fully  stressed 
design;  but,  nevertheless,,  the  results  indicate  the  apparent 
effectiveness  of  the  optimization  and  analysis  combination. 


Revolved  Section 

I 


P  =  3000  .lb 
(Distrib. ) 


Figure  2.  Two-Dimensional  ‘‘Beam’*  Problem 


TABLE  1 


ISOTROPIC  BEAM-PLATE  OPTIMIZATION 


Design 

Variable 

Initial  Design 

Final  Design 

Y1 

5.000  in. 

1.845  in. 

y2 

5.000 

1.055 

Y3 

5.000 

1.053 

Y4 

5.000 

0.988 

Y5 

2.500 

0.923 

Y6 

2.500 

0.527 

Y7 

2.500 

0.527 

Y8 

2.500 

0.495 

*1 

0.100 

0.0289 

t2 

0.100 

0,0289 

Following  successful  check-out  of  the  optimization  and  basic 
analysis  logic,  the  analysis  routine  was  expanded  to  include  all 
logic  associated  with  composite  structures.  Included  were  the 
computation  of  equivalent  homogeneous  composite  properties  from 
the  basic  layer  (the  stacking  process)  end  a  procedure  for 
transforming  gross  material  stresses  into  individual  layer 
stresses  for  determination  of  failure  modes.  In  addition,  and 
under  company  support,  a  rectangular  anisotropic  finite  element 
was  developed  for  use  in  the  overall  analysis  procedure. 

check-out  was  initiated,  but,  due  to  unexpected 
difficulties  and  computational  expense  associated  with  the  finite 
element  analysis  routines,  verification  was  not  completed,  in 
spite  of  having  introduced  a  simplified  analysis  approach.  No 
production  executions  of  the  program  were  therefore  accomplished 
The  following  paragraphs  describe  effort  that  was  completed  cm 
the  program. 
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Stress  Analysis  Procedure 


Finite  element  methods  of  structural  analysis  are  widely 
used  because  of  their  advantages  in  generality  of  application  to 
structural  forms.  With  a  large  library  of  element  types 

developed  for  the  analysis  of  many  types  of  structures,  analysis 
can  be  readily  accomplished  in  a  logical  and  automated,  manner. 
The  representation  of  the  full  equivalent  homogeneous  anisotropic 
material  properties  for  assembled  multilayer  fiber  composites 
(including  effects  of  membrane-flexural  coupling)  is  carried  out 
using  a  stacking  process.  Each  layer  is  modeled  by  an 

orthotropic  parallelepiped  with  principal  axes  corresponding  to 
he  fiber  direction.  By  application  of  the  condition  of 
Linearity  of  a  surface  normal,  the  resultant  multidegree-of- 
freedom  model  is  condensed  to  yield  a  stiffness  tensor  for  the 
assembled  composite  material  in  terms  of  gross  membrane  forces 
;nd  moments . 

For  the  design  optimization,  a  flat  panel  stiffened  by  hat- 
section  stringers  was  selected  (Figure  4  )  .  The  basic  repeating 
section  of  the  panel,  after  application  of  symmetry  conditions, 
was  modeled  for  a  finite  element  analysis  by  18  rectangular 
nlates  as  indicated  in  Figure  5.  In  order  to  include  full 

membrane-flexural  coupling  effects ,  a  new  rectangular  finite 
element  was  developed  for  use  in  Bell's  General  Purpose  Analysis 
Program.  The  element  was  a  simplification  of  a  cylindrical 
rectangular  element  derived  initially  in  Reference  15. 
Formulation  of  the  element  relationships,  stiffness,  and  stress 
matrices  was  very  complex,  even  with  elimination  of  terms  in  the 
matrices  associated  with  cylindrical  curvature.  The  element  was 
successfully  coded,  and  check-out  procedures  were  initiated.  At 
this  stage,  it  was  found  that  excessively  long  computer  times 
were  required,  Foi  the  simple  generation  and  assembly  of  a 
single  element  stiffness  matrix,  not  including  the  stress  matrix, 
approximately  2  minutes  of  IBM  360/65  computer  time  were 
required.  For  the  full  18-element  model  for  the  stiffened  plate, 
the  assembly  time  was  over  30  minutes.  Solution  of  the  resultant 
equations,  assembly  of  stress  matrices  and  final  decomposition  of 
stresses  down  to  layer  level  increased  the  total  time  foi  a 
single  analysis  to  over  1  hour  on  the  360/65  computer.  Even  for 
analysis  these  times  and  associated  r  sts  are  excessive,  while 
for  an  optimization  program  involving  many  analyses  during  the 
search  procedure,  the  program  becomes  economically  unfeasible. 

From  the  past  operational  experience  with  computer-based 
analyses,  the  analysis  of  a  large  complex  structure  (e.g.  15,000 
degrees-of-freedom)  involves  approximately  30  hours  of  computer 
time.  This  condition  may  be  regarded  as  a  practical  upper  limit 
on  a  single  analysis  with  current  computer  capabilities.  Larger 
running  times  would  bring  important  questions  of  computer 
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Figure  4. 


Original  Flat  Panel  Section 


Figure  5. 


Stiffened  Panel  Geometry  and  Idealization 


reliability  into  major  consideration.  The  use  of  a  finite 
element  analysis  model  for  optimization  of  complex  composite 
structures  which  require  over  1  hour  per  cycle  (where  hundreds 
of  cycles  are  involved)  must  therefore  be  regarded  as 
economically  unacceptable .  The  technical  feasibility,  however, 
is  in  no  way  challenged  by  this  consideration. 

A  number  of  attempts  were  made  to  effect  some  degree  of 
simplification  in  the  analytical  formulation  of  the  element 
matrices,  but  no  significant  reduction  could  be  accomplished.  It 
was  therefore  considered  necessary  to  discard  the  finite  element 
approach  and  to  substitute  a  less  rigorous  analytical  approach  to 
increase  computational  efficiency. 

A  modified  analytical  approach  was  therefore  selected  which 
could  be  applied  to  a  representative  itodel  of  the  fuselage 

component.  For  the  modified  analysis,  engineering  bending  theory 
was  used.  The  structural  model  was  taken  as  a  full  circular 
cylinder  stiffened  by  equally  spaced  hat-section  stringers 

(Figure  6  ) .  It  was  assumed  that  one  stringer  is  situated  at  the 
top  dead  center  of  the  cross  section.  All  stringers  were  taken 
as  being  identical  and  fabricated  from  the  same  material  layup. 
The  outer  skin  layup  could  be  different  from  the  stringer  layup. 

The  cylinder  was  analyzed  for  three  different  types  of 
applied  loads:  bending  moment  about  a  horizontal  axis,  vertical 
shear  force,  and  torque.  The  moment  gives  rise  to  direct 

membrane  s  ^essas,  while  shear  and  torque  each  produce  shear 
stresses  in  skin  and  stringers.  In  order  to  simplify  the 
analysis,  several  basic  assumptions  were  made  as  follows: 

(a)  There  are  many  stringers  and  their  cross-sectional 

dimensions  are  small  compared  with  the  cylinder 
diameter. 

(b)  Boundary  effects  are  ignored, 

(c)  The  skin  is  fully  effective  in  both  bending  and  shear. 

Since  all  stringers  have  identical  sections,  the  complete 
stress  distributions  were  computed  for  one  stringer  only.  These 
are  given  in  Table  XI.  The  dimensions  of  the  stringer  are 

given  in  Figure  7.  There  are  n  stringers  in  the  complete 

cylinder  and  all  stresses  are  quoted  for  the  ^-th  stringer.  The 
stringer  at  the  top  dead  center  is  stringer  0  but  must  be  treated 
as  the  n-th  when  calculating  stresses  according  to  Table  II. 

The  calculation  of  stresses  using  this  procedure  is  extremely 
rapid.  The  decomposition  process,  which  is  a  reverse  of  the 
stacking  procedure,  is  still  applied  after  calculation  of  gross 
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Figure  6.  Cylinder  Cross  Section 
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Figure  7.  Detail  of  Single  Skin  -  Stringer 
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stresses  to  determine  individual  layer  stresses  for  use  in  the 
failure  criteria  expression. 

The  overall  program  logic  presented  in  Reference  3  remains 
substantially  unaltered  by  this  new  formulation.  The  program 
modules  associated  with  analysis  of  the  finite  element  model  have 
been  removed  completely  and  replaced  by  a  simple  routine  which 
evaluates  the  overall  stresses  according  to  the  bending  theory 
analysis  of  Table  II.  All  logic  prior  to  and  following  the 
analysis  was  effectively  unaltered  by  this  substitution. 
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SECTION  IV 


BEHAVIOR  OF  COMPOSITE  STRUCTURAL  ELEMENTS 


A.  Post-Buckling  of  Anisotropic  Plates  Under  Axial  Compression 
(Professor  Richer  and  Mr.  D.  Chan) 

The  post-buckling  behavior  of  laminated,  anisotropic  plates  was  inves¬ 
tigated.  The  Love-Kirchoff  assumptions  for  thin  plates  and  Von-Karman's 
large  displacement  kinematic  relations  were  used.  The  potential  energy  ex¬ 
pression  was  formulated  on  the  basis  of  an  extensional  plate  theory.  Plates 
treated  are  rectangular  with  all  four  edges  either  clamped  or  simply  supported. 
The  plates  are  loaded  on  two  opposite  sides  with  the  unloaded  edges  either 
laterally  restrained  or  free  to  expand.  The  loaded  edges  are  subjected  to 
either  uniform  force  or  uniform  displacement.  All  load  levels  are  such  that 
the  plate  material  remains  elastic.  The  solution  was  carried  out  along  the 
line  of  Raleigh-Ritz  method  with  only  the  transverse  displacement  mode,  w, 
assumed;  the  remaining  two  displacement  modes,  u  and  v,  being  obtained  .exactly 
from  the  two  nonlinear  membrane  equilibrium  equations. 

The  results  indicate  that  uniform  force  load  and  uniform  displacement  load 
yield  exactly  the  same  post-bu kl ing  load  vs.  deflection  curve.  For  most  plates, 
buckling  is  a  bifurcation  phenomenon;  however,  for  certain  types  of  plates  with 
certain  types  of  boundary  conditions,  tho  plate  deflects  as  soon  as  the  load  is 
applied.  Solutions  were  obtained  for  simply  supported  plates  with  a  single 
transverse  mode  assumption  and  with  a  double  transverse  mode  assumption. 

* 

The  method  of  solution  has  been  put  into  a  Fortran  program  which  requires 
the  plate  aspect  ratio,  thickness,  material  constants  and  wave  numbers  as  in¬ 
puts  and  yields  the  load  vs.  deflection,  shortening,  effective  width  and  re¬ 
sultant  force  distribution  relations  as  outputs.  Average  computing  time  on 
Uni  vac  1108  for  plate  analysis  is  between  four  to  seven  seconds. 

Figures  8,  9,  and  ’0  present  the  load  deflection  curve,  effecting  width 
versus  shortening  curve  and  the  load  shortening  curve  for  a  (0,90,90,0) 

'‘Thornel"  -  25  plate  respectively.  This  plate  was  tested  earlier  in  the  pro¬ 
gram  by  Mandell  (3)  and  analyzed  by  Monforton.  Figure  8  shows  that  the  results 
compare  favorably  during  the  early  stages  of  loading.  Improvements  in  the 
post-buckling  prediction  could  be  achieved  by  assuming  more  terms  for  the  w 
displacement  expression. 

A  topical  report  on  this  subject  has  been  completed  and  submitted  to  Case 
Western  Reserve  University  by  D  P.  Chan  as  partial  completion  of  the  degree 
requirements. 


B  Post- Buckl ing  Strength  of  Composite  Plates 

(Professor  T  P  Kicher  and  Mr.  M.  T.  Islam) 

Recent  studies  on  the  buckling  and  (>ost  buckling  of  laminated  plates  in¬ 
dicate  that  significant  loads  can  be  sustained  after  buckling  if  the  associated 


-  22  - 


Figure  8.  Load-Deflection  Curve  for  (0,90,90,0)  "Thornel"  25  Plate 
Simply  Supported  and  with  UnLoaded  Edges  Free  to  Expand 
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Figure  9.  Effective  Width  vs.  Shortening  Curve  (0,90,90,0)  "Thorne!"  25  Plate 
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clamped  plate 


Figure  10. 


Loao  Shortening  Curve  for  (0,90,90,0) 


"Thornel "  25  Plate 
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larqe  displacements  can  be  tolerated.  The  high  stiffness  of  composites  leads 
to  thin  laminates  for  lightly  loaded  structures,  consequently, the  high  strength 
potential  is  not  utilized.  If  analytical  methods  of  predicting  the  post-buckling 
response  and  strength  were  available,  designers  could  utilize  the  full  potential 
of  the  post-buckling  region.  The  post-buckling  response  is  reported  in  the 
previous  section  while  the  strength  analysis  is  reported  herein. 

1 •  General  Edge  Loading 

The  strain  displacement  relations  for  plates  under  large  deflections 
can  be  given  by 

o 

£x  =  ex  +  ZKX 


ey  +  ZKy 


(l-i) 


Y  =  Y  +  Z  K 

Yxy  'xy  xy 

where  Kirchoff-Love  conditions  are  assumed  and  z  is  measured  from  the 
reference  surface.  The  reference  surface  strains  are  given  by 

o  ]  2 

ex  =  ux  +  2*x 
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Yxy  uy  +  ux  +  wx  wy 
and  the  curvature  terms  are 
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The  subscripts  of  u,  v  and  w  terms  indicate  partial  derivatives.  The  force 
deformation  relationship  for  the  general  anisotropic  plate  is  then 
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(1-6) 
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where  s  indicates  the  surface  of  the  plate.  Then  the  strain  energy  of  the 
plate  is 


U  =  U]  +  u2  +  u3 


(1-8) 


The  potential  of  the  applied  loads  is 


W  = 


[  "x  “  4  "x,  v  *  ■  \  wx  -  Bxy  wy]  d* 
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where  Nx,  Rxy,  Nyx,  Qxz>  Oyz  are  forces  and  Ry,  Ry,  Rxy,  Ryx  are  moments 
applied  on  the  edges  of  the  plate. 

The  total  potential  energy  np  of  the  plate  is:  np  =  U  -  W  (1-10) 

For  prescribed  kinematical  boundary  conditions,  states  of  equilibrium  of 
the  plate  will  correspond  to  stationary  values  of  the  elastic  energy.  Under 
hinged  or  built-in  edge  conditions  the  elastic  energy  consists  of  the  energy 
of  the  plate  only.  For  elastically  restrained  plate  edge  conditions,  the 
elastic  energy  of  the  edge  stringers  must  be  added  to  obtain  the  total  elastic 
energy.  In  any  case,  the  stationary  value  of  the  elastic  energy  may  be  cal¬ 
culated  s  a  function  of  the  applied  deformation  l.  The  approximate  solution 
may  be  obtained  by  Rayleigh's  or  Rayleigh-Ritz's  principle.  Let  the  solution 
obtained  by  u(e)»  V(s),  w(e)  for  an  applied  deformation  l  .  Call  the  elastic 
energy,  which  is  a  function  of  l  ,  U(e). 

Then 
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0-11) 


dU  = 
or  —  =  F 


where  F  is  the  edge  load  (vector).  The  ultimate  load  capacity  is 


e  _  max  dU 
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2.  End  Compressive  Loading 

,.nWD(,  a  „  Consid<:r  an  orthotropic  or  anisotropic  plate  of  length  a  and  width  b 

^  L a  c?mPre?s1ve  1oad  p*  The  load  is  applied  through  a  ri;id  bar  at  the  end 
and  the  plate  is  supported  on  all  sides.  ’ 

n^cc?Urin?  t-6  fundamental  or  unbuckled  state,  the  plate  suffers  uniform  com- 

alonu X  w1d‘h-the”‘^  l°*d  is  distributed  uniformly 

along  the  edge.  After  buckling,  the  center  of  the  plate  bows  out  while  the 
siues  remain  straight.  Because  of  the  membrane  stretching  caused  by  the  de- 
-  r  at  t(|e  middle,  the  compressive  stress  there  is  reduced.  To  support 
the  same  applied  load,  the  sides  of  the  plate  then  carry  proportionately  greater 
t0*?1  I"  «»  far  post-buckled  state,  the  load  JistriSlon 

wLth  h°Wn  1n  FlguJ"e  ]J^a)  where  the  entire  load  is  carried  by  two  strips  of 
width  d,  one  on  each  side.  The  loading  is  considered  to  be  uniformly  distri¬ 
buted  across  these  two  strips  each  carrying  half  the  total  load  P.  Let  the 

momKrano  rnmnrorrl,,.  1 _ I  __  J_| _  ,  .  .  ..  ’  r  *  t-CL 


membrane  compressive  load  on  these  strips  be  N  .  Then 


P  =  2dN 


(2-1) 


I*®  plate ,1S  assumed  to  be  simply  supported  on  all  sides.  The  lateral 
deflectTon  w  along  the  edges  is  thus  zero.  For  the  purpose  of  analysis  the 
middle  strip  without  edge  leading  is  ignored.  (Fig.  12)  y  *  tne 

Orthotropic  Case 

trnmV1  nif!cal . ?['thoJroP1  ^  buckling  theory  can  now  be  applied  to  the  ortho- 

Jf?  with  material  axes  aligned  with  structural  axes  of  the  plate,  and 
with  the  plies  oriented  synmetrical ly  about  the  mid-plane  of  the  plate.  This 
b^eory  gives  the  lateral  deflection  w  due  to  compressive  load  along  the  x  axis 
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where  D..  are  flexural  stiffnesses  and  D,  =  + 
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A  solution  of  equation  (2-2)  is  given  by 


W  =  w  sin  ~  sin 


(2-3) 
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Figure  11 

Post-buckling  Load 
Distribution 


Figure  12 

Post-buckling  Analysis 
Model 
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where  a  is  the  length  of  a  half  wave  along  the  x  direction.  Equation  (2-3) 
satisfies  the  boundary  conditions  of  the  plate  i.e.,  w  =  0  and  w  =0  at 

r  =  0  or  A  and  w  =  0  and  w>yy  =  0  at  y  =  0  or  2d.  Substituting  (2-3)  into 

(2-2)  and  dividing  through  by  W  sin  ~  sin  ^  we  get 
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We  assume  the  two  strips  carrying  the  load  are  at  the  stability  limit 
corresponding  to  the  critical  wave  length  a  .  Differentiating  the  terms  inside 
the  square  brackets  with  respect  to  a  and  equating  it  to  zero 


Substituting  in  (2-4) 
2 


N 


x  [  V  D11  °22 


d  *  jir  I  D"D22 


and 


P.2dNx=^*  C(/D 


11 


Replacing  Nx  by  at 


d  = 


V 


22 


p  =  ft  *  [(  /d^d 


22 


+  03  3 

(2-5) 

1  1/2 

+  °3J 

(2-6) 

°22  +  °3*  Nx^/2 

(2-7) 

-1  1/2 

+  D3j 

(2-8) 

+  D3)  oxt],/2 

(2-9) 

Equations  (2-8  and  2-9)  give  the  relation  between  effective  width  and  the 
load.  These  equations  hold  up  to  the  elastic  limit  or  yield  point  for  the 
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material.  However,  they  may  be  used  beyond  the  elastic  limit  if  D^,  D22  and 

Dg  are  computed  on  the  basis  of  the  tangent  moduli  for  the  material  of  the  plate. 

For  composites  where  the  material  if,  elastic  up  to  the  stress  corresponding  to 
failure  the  ultimate  load  capacity  Is 

P„lt  -V*-'  [<  'TTO  +  D3)»xc«,/2  (2-,0) 

where  0  is  the  failure  stress  in  longitudinal  compression  for  the  composite, 
xc 

Unbalanced  Case 

In  the  case  of  an  orthotropic  plate  whose  laminae  are  not  elastically 
symmetric  with  respect  to  the  reference  surface,  the  plate  will  have  bending- 
membrane  coupling.  The  governing  equation  for  deflection  for  such  a  plate  is 


M  +2M  +  M  +  N  W  =0 

a  ,xx  xy,xy  y,yy  x  ,xx 


(2-11) 


along  with  the  force  deformation  relationship 


|{M} 


>  = 


[A] 

[B] 


[BJ 

CD] 


(el 

M 


(2-12) 


The  solution  for  {M}  from  equation  (2-12)  is 


{M>  =  [D  -  BA-1B]  {<}  -  [BA-1]  {N}. 


(2-13) 


For  simplicity  neglecting  the  term  [BA'1]  {N>  in  equation  (2-13)  leads 
to  the  governing  equation 

D11  w’xxxx  +  2D3  W’xxyy  +  °22  W,yyyy  +  Nx  W,xx  0  ^  14^ 

where  the  D*  terms  are  elements  of  the  matrix  [D-BA"1  B]. 

Equation  (2-14)  is  of  the  same  form  as  (2-2),  and  analogus  treatment  as  in 
the  orthotropic  case  gives 


d  = 


/ 


*  *  * 
D11  °22  +  °3 


1/2 


(2-15) 
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=  rr  * 


[(  / D*u  0*2  +  D*)  oxt  j1 


(2-16) 


*h<,Uit1!!ns  f2"-b  2'16^  9ive  the  relatlon  between  the  effective  width 

nd  the  1oJd-  Js  in  the  orthotropic  case  these  equations  hold  up  to  the  elastic 

limit  if  D^,  022  and  D3  are  computed  on  the  basis  of  tanqent  moduli.  For  com- 
?ap^trrialS  WhiCh  remain  e1aStiC  UP  t0  tailure’  we  obtain  the  ultimate  load 


=  rr  n  [(  /d*u  022  +  D*  )  a  t] 


(2-17) 


where  axc  is  the  failure  stress  for  the  anisotropic  composite  in  longitudinal 
compression. 

Rations  (2-10  and  2-17)  have  been  used  to  compute  the  ultimate  load 
capacities  of  various  composite  plates  with  different  ply  lay-ups  and  thick- 

i^5eS'n  Rl-  tS-  indicate  considerable  reserve  strength  beyond  the  bucklinq 
load.  Work  is  in  progress  towards  verification  of  the  theories  throuqh 
experiment.  ^ 


Buckling  of  Laminated  Composite  Plates  with  Cutouts 
(Professor  T.  P.  Kicher  and  Mr.  J.  Martin) - 


^The  purpose  of  this  study  is  to  investigate  the  buckling  of  laminated  com- 
posite  plates  with  cutouts,  and  to  determine  the  effect  of  the  cutouts  on  the 
duck i ing  i odds . 

in  ,’n"  rlrna7  was  8-ply  laminated  square  plate,  . 

in  tho  in*Tuy  ‘?8?  in”  W1,h  2  in<  dnd  4  ln-  diameter  circular  cutouts 

inn9heaCen!'3rui  The  plat?  was  mflde  from  "Scotchply"  Reinforced  Plastic  Type 
on  on**  m?cda5Ie  ?°r9  ??s  larmnate-  The  ply  orientations  were  0,  45,  -45, 
.'V  45»  5*  Buckling  loads  were  found  experimentally.  The  plate  was 

Hn?  ne  Jho1  i  simP  y-supported  on  all  four  edges.  Microscopic  inspec¬ 
tion  of  the  plate  revealed  large  and  numerous  voids. 

HowoiJ"  3rder  t0  Produce  uniform  laminated  plates  a  fabrication  process^/as 
developed  using  a  Patterson-Kelly  autoclave.  Two  10-ply  unidirectional 

"Sco?rhnt'qpfrnfP  !S^i12j‘n-rby  12  in*  b^  -100  in-»  were  ^Heated  from 
wr.v/D^Chf  y  Reln!°l?ed  Plastlc  Type  1002.  A  vacuum  bag  technique  was  employed- 
woven  glass  and  Teflon  were  used  between  the  uncured  laminate  and  a  piece  of 

Tf  p,?ce  on  top  of.the  fi™'  'w  Of  woven  fllass.  Gene  nresoS-e 

It  330  ntai"ed  at  18-?°  r,s1>  and  texture  «as  maintained 

rnmniot  i  L/1}9  the  35  ™nute  cure  cycle.  After  the  cure  cycle  had  been 

d>  the  p1ate,was  allowed  to  cool  very  slowly  inside  the  autoclave  while 
tne  gage  pressure  and  vacuum  were  still  maintained.  This  procedure  produced 
a  plate  wmch  was  very  flat.  Fiber  wash  was  found  to  be  negligible,  and  void 
content  was  low.  The  average  cured  thickness  was  .010  in.  per  ply. 
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brittle  coating  techniques  were  applied  to  a  laminated  plate  havinq  a 
circular  cutout  in  a  preliminary  investigation  of  the  strain  field  for  the 
plate  under  axial  load. 

A  literature  search  was  started  pertaining  to  buckling  of  laminated  com¬ 
posite  Dlates  with  cutouts;  it  saems  as  though  the;  published  literature  on  this 
topic  is  very  small. 


D ■  Accuracy  Study  of  Shell  Theories  for  Anisotropic  Cylinders 

(Professor  T.  P.~  Kicher  and  Mr.  C.  !l.  Wu) 


The  major  work  of  the  investigation  of  shell  theories  for  anisotropic  cy¬ 
linders  was  completed  during  this  report  period.  The  following  are  brief 
descriptions  of  the  work  done: 


(a)  Formulation  of  boundary  determinants  for  buckling:  Three  sets  of  dif¬ 
ferential  equations  are  used  for  the  analysis  of  shell  buckling.  They  are 
based  on  the  theories  of  Fluggeu  Donnell,  and  Timoshenko.  Eight  boundary  con¬ 
ditions  are  const  cere'';  four  simple  support  cases  and  four  clamped  end  cases. 
For  example: 


Simple  support  SS  I: 


Clamped  End  Cl: 


w  =  0 
Hx  -  0 

v  t  0;  or  Tx  +  ?  v,x  *  0 
u  f  0;  or  Nx  +  j-  u>(J  =  0 


w  =  0 
w,x  =  ° 

v  i<  0;  or  Tx  +  P  v,x  »  0 
u  f  0;  or  Nx  +  j  u,0  *  0 


Mx 

where  Ty  =  Ny  +  —  is  the  Kirchoff  boundary  force.  The  other  three  cases  can 
be  obtained  by  setting  u  =  0,  v  =  0  or  u  =  0  and  v  s  0. 


The  following  displacement  modes  are  assumed: 

8  A.  x 

u  =  z  U.  sin  (  +  ne  ) 

k=l  k  a 


v 


A.X 

sin  (  +  ne  ) 

a 


w 


w  A.  A 

z  W.  cos  (  -~  +  ne  ) 
k=l  k  3 


(1) 
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Wh<ire  V  Uk’  V  and  wk  are*  in  general,  complex  values.  The  assume  HicP?ace 

the  boundarv^nnrfit •4atl°nS  but  the  bounda^  conditions, 
formulated.  Sis  is  done  by  express fna°thp  fnif6C-al  bounda7  determinant  is 
imaginary  parts:  y  P  eSS1"9  the  Allowing  quant  ties  in  real  and 

Ak  =  (Ar  +  iAi  )k 
Wk  =  <Wr  +  «i  )k 

Uk  *  <Ur  +  iUi  )k  -  <  f„r  "  tfut  >k  Wk  ,, 


Vk  =  <Vr  +  \ 


^  fvr  +  lfvi  ^  wk 


(2) 

quanti ^respectively!'  ^  1  den0te  the  real  and  the  ima9inary  P^t  of  the 

Eqs.  Eq;-  0). 

Into  the  appropriate  bou£dar.TSndU*ons  Fo?  w  -  n  f?rce-stra’n  relations 
pressions  for  the  boundary  det^Ln?  ^  "'  °  the  C°"-«P°"‘»"9  ex- 

7 

k_7.3  (“kWsrk  *  Vl  “clktl  >  ■  ° 

7 

k=l,3  ^kWsrk  +  “k+1  Wsik+1  )  =  0  if  a  complex. 


or 


k 


8 

l 

k=l 

8 

l 

k=l 


ukwcrk  =  0 


ukWsrk  =  0 


if  x  real. 


S1milar  expressions  can  be  obtained  for  u-0,  v=0,  w,  =0,  M  =0,  T  +  P  v  -  0 

X  X  x  *x  ~  u’ 

and  Nx  ♦  I  u.e  -  o. 

linear°equat1ons  2aVb1  ?’?»  „ 

For  nontnml  solution  of  v  the  boundary  determinant  must  equal  aero.'  The®' 

cylinder  Wave  "•  *r  a 

which  satisfy  the  vanishing  of  the  Lnda™  de'ie^n?"*5  those  va,ues 
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(b)  Computer  programs:  The  computer  programs  for  all  three  theories  of 
shell  buckling  have  been  completed.  After  an  extensive  testing  and  study  of 
the  programs  and  the  numerical  results,  it  is  found  that  the  iteration  proce¬ 
dures  give  very  good  results  for  the  case';  of  radial  pressure,  torsion  or  a 
combination  of  both.  Figure  13  shows_the  interaction  curve  for  the  combined 
radial  pressure  and  torsion  with  4(0,+  22.5,  90)  layer  configuration,  SS4 
boundary  condition,  L/R=2  and  R/t=l 61 .6.  FlUgge's  equations  were  used  in  this 
case.  A  similar  curve  wis  obtained  for  Donnell's  theory.  The  differences  be¬ 
tween  Fliigge  and  Donnel  i  are,  in  general,  less  than  1%.  The  interaction  curve 
resembles  the  shape  of  a  parabola.  It  is  observed  that  the  internal  pressure 
tends  to  strengthen  t  cylinder  while  the  external  pressure  tends_to  weaken 
the  cylinder.  Because  of  the  unbalanced  layer  configuration,  4(0, +22. 5,  90), 
the  curve  is  slightly  deviated  from  symmetry  about  the  horizontal  axis. 

The  results  for  axial  compression  obtained  by  the  usual  iteration  proce¬ 
dures  are  merely  the  approximate  solutions.  The  difficulty  in  obtaining  an 
exact  buckling  load  for  the  axial  compression  case  comes  from  the  basic 
characteristic  of  the  critical  loads.  The  critical  loads,  or  the  eigenvalues, 
for  the  axial  compression  case  lie  very  closely  togetner.  By  using  the  usual 
iteration  procedures  with  coarse  load  increments,  the  computer  program  might 
miss  the  lowest  eigenvalue  and  converge  to  a  higher  one.  To  overcome  this 
difficulty  it  is  necessary  to  use  the  approximate  buckling  load  and  a  very 
f!ne  load  increment  to  pinpoint  the  lowest  buckling  load.  This  modified 
ueration  procedure  has  been  used  with  considerable  degree  of  success.  One 
major  drawback  of  the  procedure  is  the  enormous  amount  of  computer  time 
consumed. 

(c)  Static  responses:  Equations  and  methods  of  solution  for  the  static 
responses  due  to  the  surface  loading  were  formulated.  Computer  programs  have 
been  completed.  The  numerical  results  are  currently  under  study. 
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FIGURE  13.  Radial  Pressure  and  Torsion  Interaction  Curve 
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SECTION  V 


MICROMECHANICS  OF  FIBER  COMPOSITES 


A .  Three=Dimensional  Photoelastic  Models  of  Filamentary  Composite  Materials 
(Professor  D.  K.  Wright  and  Mr.  D.  Weitzenhof) 

The  development  of  a  suitable  photoelastic  technique  for  analysis  of  fila¬ 
mentary  composites  has  been  completed.  Tests  show  that  three  dimensional  stress 
states  can  be  evaluated  with  reasonable  accuracy  using  photoelastic  models. 

The  frozen  stress  method  of  loading  was  selected  in  conjunction  with  a 
strain  energy  failure  theory.  This  approach  avoids  errors  introduced  by  in¬ 
tegration  techniques  such  as  the  shear-difference  method. 

1 .  Development 

For  east  in  analysis  and  handling,  the  frozen  stress  method  of  photoelastic 
analysis  was  selected.  For  simplicity,  all  models  were  machined  and  cemented 
together  using  PLM-1  clear  model  cement  from  Photolastic,  Inc.  The  most  formid¬ 
able  obstacle  to  the  accurate  use  of  this  procedure  was  the  formation  of  large 
thermal  stresses  which  were  usually  of  the  same  magnitude  as  the  stresses 
inuuced  by  external  loading.  These  large  thermal  stresses  required  the  use  of 
two  models,  a  loaded  model  and  an  externally  nonloaded  model.  Each  model  had 
to  be  analyzed  for  the  complete  stress  state  at  each  point  in  question.  These 
stress  states  were  then  subtracted  in  accordance  with  superposition  theory  to 
give  the  stress  state  produced  by  the  external  loading.  Since  the  completa 
stress  state  had  to  be  known  and  subtracted,  t.he  direct  use  of  strain  energy 
theory  or  similar  approaches  was  impossible. 

Since  these  thermal  'stresses  were  caused  by  the  difference  in  thermal 
expansion  coefficients  between  the  photoelastic  matrix  model  or  the  metal  or 
wooo  fiber  model,  the  obvious  solution  was  to  use  some  fiber  model  material 
with  a  more  compatible  coefficient  of  thermal  expansion.  To  this  end,  and  to 
reduce  the  magnitude  of  the  thermal  stresses  still  further,  an  extremely  low 
critical  temperature  (170°F)  photoelastic  epoxy  was  developed  by  Photolastic, 
Inc.  This  low  critical  temperature  allowed  the  use  of  a  plastic  for  the  fiber 
material.  Two  plastics  selected  for  testing  were  Plexiglas*  and  Lexan.** 

Preliminary  tests  were  conducted  on  the  photoelastic  epoxy,  Plexiglas,  and 
the  Lexan  to  determine  their  properties  at  170°F.  Tests  were  also  conducted 
on  Lexan-epoxy  and  Plexiglas-epoxy  models  to  check  for  fiber-matrix  compati¬ 
bility  during  stress-freezing.  The  material  properties  tests  gave  results 
as  shown  in  Table  III. 


*  Rohm  &  Haas  Co. 

**  General  Electric  Co. 


TABLE  III 

PROPERTIES  OF  PHOTOELASTIC  MATRIX  MATERIALS 


MATERIAL  ~ 

VOUNG'S 

MODULUS 

"  POISSON'S 
RATIO 

»  COMMLNIS 

Photoelastic 

Epoxy 

1890  PS I 

.48 

Full  Deformation  Reached  In  Six 
Minutes 

Lexan 

250  KSI 

.35 

Modulus  Measured  After  24  Hrs 

At  170°  F 

Plexiglas 

5  to  150 

KSI 

.40-. 44 

Excessive  Creep  at  170°F 

Tabled  md,lS  (,Fi9UreS  14_16)  Produced  results  as  shown  in 

om  these  results,  Lexan  was  selected  as  the  fiber  model  material. 


Tests  7  and  8  were  then  conducted 
the  thermal  and  load- induced  stresses 
fringe  patterns  produced.  These  tests 
method  discussed. 


to  determine  the  relative  values  of 
Figures  17-21  show  configurations  and 
demonstrated  the  feasibility  of  the 


TABLE  IV 


TEST  RESULTS  ON  VARIOUS  COMPATIBILITY  MODELS 


Test 

NO. 

MODEL 

TYPE 

“Loading  ' — 

“Maximum 

SHEAR  STRESS 

COMMENTS- - 

1 

2-D 

LEXAN 

THERMAL 

ONLY 

2.9  PSI 

FOR  CONFIGURATION 

SEE  FIGURE  14 

2 

2-D 

PLEXIGLAS 

thermal 

ONLY 

1.0  PSI 

SEE  FIGURE  14 

3 

2-D 

LEXAN 

2.63#  TENSILE 
PLUS  THERMAL 

8.5  PSI 

SEE  FIGURE  15 

4 

2-D 

PLEXIGLAS 

2.63#  TENSILE 
PLUS  THERMAL 

8.1  PSI 

SEE  FIGURE  15 

5 

3-D 

LEXAN 

THERMAL 

ONLY 

5.4  PSI 

SEE  FIGURE  16 

6 

3-D 

PLEXIGLAS 

THERMAL 

ONLY 

5.4  PSI 

SEE  FIGURE  16 

7 

3-D 

LEXAN 

THERMAL 

ONLY 

12.0  PSI 

FOR  CONFIGURATION  SEE 
FIGURE  17.  FOR  FRINGE 
PATTERNS  SEE  FIG  is  *  iq 

8 

3-D 

PLEXIGLAS 

152#  COMPRESSIVE 
PLUS  THERMAL 

58.8  PSI 

FOR  CONFIGURATION  SEE 
FIGURE  17.  FOR  FRINGE 
PATTERNS  SEE  FIG. 20. &  21 

-  39  - 


-  40  - 


FRACTION  ♦ 


Figure  15. 


Model  Type  2  - 


1"  Thick  Full  Scale 

J 
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Low  Temperature  Epoxy,  Lexan 
and  Plexiglas  Fibers 


Section  A  -  4 


Figure  16.  Preliminary  S.  F.  Test  Model  Sketches 
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1" 

T 


Section  A  -  A 


Figure  17.  Load  Test  No.  1  Model  Sketches 


Low  Temperature  Epoxy, 
Lexan  Fibers 
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Figure  18.  Thermal  Stresses  Near  a 
Fiber  End  in  a  Dark  Field. 

(1.25  Fringe) 

N-23037 
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Figure  19.  Thermal  Stresses  Near 
Fiber  End  in  a  Light  Field. 
(1.25  Fringes) 


a 

N-23038 
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Figure  20.  Thermal  and  Mechanical  Stresses  Near  a 
Fiber  End  in  a  Dark  Field. 

(6.0  Fringes)  N-23039 
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Figure  21  . 


Thermal  and  Mechanical  Stress 
Fiber  End  in  a  Light  Field. 
(6.0  Fringes) 


es  Near 
N-23040 


a 
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2.  Strain  Energy  Analysis 


The  strain  energy  equations  may  be  stated  as  follows  (15): 

where 

tAB  tAC 

°B  tBC 

tCB  aC 

and  E  is  the  equivalent  strain  matrix.  Eli mi  rating  the  hydrostatic  stress 
component  which  does  not  contribute  to  failure,  and  realizing  that 


°A 

tBA 

tCA 


E 


S 

E 


where  E  is  Young's  Modulus  we  obtain 


.  (°r°2)2  *  (o2'°3)2  *  (°3~°1)2 
126 

where  G  =  20+zrr  • 

9  2  2 
2o,  o 

In  simple  tension,  U'  =  =  gf-  • 

Thus,  failure  will  occur  when 

(a- |-a2)2  +  (°2_a3^2  +  ^a3"cl^  =  ^ao 
Using  photoelastic  relationships,  this  reduces  to 


2 

aeff 


l  (^)2  0  +  jsin2  2<f>) 


where  f  is  the  stress-fringe  constant  for  the  particular  material;  t  is  the 
size  of  the  cube;  N  is  the  fringe  number  in  that  view  of  the  cube;  and  <p  is 
the  isoclinic  angle  in  that  view  of  the  cube. 
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2 

Since  (sin  2<f>)  is  the  same  for  (+<f>  +90°),  the  orientation  of  the  cube 
need  not  be  known.  The  only  requirement  Ts  that  all  three  directions  be  viewed 
and  the  position  of  the  cube  in  the  model  be  known. 

3.  Technique  Application 

To  demonstrate  the  technique  developed  above, several  fiber-fiber  and  fiber- 
void  interaction  models  have  been  fabricated  and  analyzed.  A  typical  model  con¬ 
figuration  and  load  application  is  shown  in  Figure  22.  Figure  23  shows  tynical 
slices  from  the  loaded  and  unloaded  models  in  the  dark  field  oolariscope.  As 
can  be  seen,  the  slice  from  the  externally  nonloaded  model  has  almost  no 
indication  of  stress  whereas  the  loaded  model  slice  shows  approximately  one  and 
one-half  fringes.  Both  slices  are  0.100  inches  thick.  The  stresses  present  in 
the  Lexan  fibers  are  caused  by  the  initial  forming  process  and  have  no  meaninq 
in  the  model.  Thus,  the  thermal  stresses,  which  were  the  greatest  obstacle  to 
a  satisfactory  analysis,  can  now  be  assumed  very  small  and  ignored.  The  cubes 
shown  were  analyzed  using  the  strain  energy  theory  discussed  above. 

4.  Results  and  Conclusions 

The  results  of  the  analysis  of  these  cubes  according  to  the  strain  enerqy 
theory  are  shown  in  Figure  24.  Other  configurations  including  voids  will  also 
be  analyzed  and  the  values  of  effective  stress  compared.  These  experimental" 
values  will  also  be  compared  to  theoretical  values  obtained  by  J.  F.  Stevenson 
in  his  work,  Discrete-Element  Microstress  Analysis  of  Unidirectional  Fihpr 
Composites.  - - — - 


Results  thus  far  indicate  that  frozen  stress  photoelasticity  may  be  used 
successfully  to  determine  with  reasonable  ease  stress  states  or  effective 
stresses  within  a  filamentary  composite  material.  A  complete  report  will  be 
forthcoming  in  the  form  of  a  thesis  submitted  in  partial  fulfillment  of  the 
requirements  for  the  degreee  of  Doctor  of  Philosophy  at  Case  Western  Reserve 
University. 

B-  Scattered  Light  Photoelasticity  and  Holography 

(Professor  J.  Cannon,  Mr.  C.  Choi,  Mr.  P.  Shankar) 

1-  Scattered  Light  Photoelastic  Studies 

The  photoelastic  studies  of  fiber  reinforced  composite  materials  to  present 
have  been  carried  out  by  the  frozen  stress  method  with  one,  two,  or  few  rein¬ 
forcing  fibers.  The  frozen  stress  method  encounters  a  serious  difficulty  of 
large  thermal  stress  induced  in  the  model  due  to  the  mismatch  of  the  thermal 
properties  of  matrix  and  fiber  materials  during  the  stress  freezing  cycle.  It 
is  believed  that  the  use  of  two  identical  models  (one  with  load  stresses  one 
without)  would  eliminate  this  difficulty. 

However,  there  still  remains  an  uncertainty  as  to  the  possibility  of 
fabricating  two  identical  models  with  multiple  reinforcinq  fibers.  One  other 
uncertainty  regarding  the  accuracy  of  the  results  with  two  models  when 
thermally  i  iduced  stress  is  much  greater  than  the  stresses  due  to  externally 
applied  load.  J 
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Configuration 


■jju 

O 

o 

■w 

Total  Load  *  250  lb. 
Load  Application 


Figure  22, 


Model  Configuration  and  Load  Application 
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•  nn?Y^j;iv®  st^y.of  th®  properties  of  transparent  plastics  show  that  it 

Dertie'  and^Hif°h  Mn  fnJ  fiber  materials  with  matching  thermal  pro- 

nf  £  w  ?d  b*  able  co  simulate  the  ^tios  of  the  mechanical  properties 
of  fiber  to  matrix  of  an  actual  composite  material.  i. 

Thus,  the  live-load  approach  to  the  photoelastic  analysis  of  the  fiber 
reinforced  composite  materials  was  adopted.  y  T  D 

Selection  of  Photoelastic  Model  Material s 
An  appropria te' combi natlo^o^the  ItrU 

forcedPcompMitetmateriaK°f  “*  ™d“H  °f  ^ 


TABuE  V 


Possible  Model 

Material 

Matrix  Material 

Thermal 

Expansion 

Coefficient 

Modulus 

Strength 

(10"5/°F) 

(105  psi) 

(10^  psi ) 

Polyesters 

Cast  Rigid 

Cast  Flexible 

3.7 

2.8  -  5.5 

1.5  -  6.5 
.001  -.10 

4  -  10 

1  -  8 

Epoxy 

Cast  Rigid 

Cast  Flexible 

3.3 

3.0  -  10.0 

4.0  -  5.0 
.01  -  3.5 

6  -  9 

2.10 

Ionomers 

6.7 

.2  -  .6 

3.5  -  5.5 

Fiber 
(or  Rod) 

Material 

Thermal 

Expansion 

Coefficient 

Modulus 

Strength 

(10"5/°F) 

( 1 05  psi ) 

(10^  psi ) 

Acryl ic 

5.0  -  9.0 

3.5  -  5.0 

8  -  11 

Polycarbonate 

3.75 

3.45 

9.5 

Aluminum 

.125 

106 

60 

An  attempt  was  made  to  determine  the  fringe  definition  characterictirc 

StfndfP^yeSter  reSi?S  °J  varying  from  rigid  to  °f 

1  aof  the  experiment  show  that  polyesters  possess  superior  fringe  deflni- 

«?!rs  SIh?Mr^SrfI°rh?Cattered-ll9ht  C°mpared  t0  epox'-  However^pofy-"1 
esters  exhibit  considerable  amounts  of  creep.  Typical  scattered-light 
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fringes  of  polyester  and  epoxy  are  shown  in  Figure  25.  The  model  is  a  disk 
uncier  diametral  compression. 

Epoxy  and  polyester  resins  used  are  as  follows: 

Epoxy 

Araldite  6020  (Resin)* 

Araldite  508  (Resin) 

Araldite  956  (Hardener) 

Polyester 

Cadco  Resin  (Rigid)# 

Laminae  EPX-126-3**  (Flexible  Resin) 

Laminae  4134**  (Flexible  Resin) 

MEK  Peroxide  (Hardener) 

Compositions  and  stress  fringe  values  are  shown  in  the  following  Tables: 

TABLE  VI 

Photoelastic  Calibration  Constants 


Sample 

f(  lb/fringe- in) 

Sample 

f(lb/fringe-in) 

Epoxy  No. 

3 

93.9 

Polyester  No.  8 

6.6 

Epoxy  No. 

4 

28.9 

Polyester  No.  10 

110.9 

Epoxy  No. 

5 

109.1 

Polyester  No.  11 

132.2 

Epoxy  No. 

6 

100.1 

Polyester  No.  12 

157.0 

Epoxy  No. 

7 

109.9 

Polyester  No.  13 

128.7 

Polyester  No.  14 

176.0 

Polyester  No.  15 

105.9 

Polyester  No.  16 

118.7 

*Ciba  Products  Company 

**American  Cynamid  Company 
#Cadi 1 1 ac  Plastic  and  Chemical  Company 
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(a)  Polyester  (#13)  (b)  Epoxy  (#7) 

Figure  25.  Scattered  Light  Fringes  for  Disks  of  Polyester 
and  Epoxy  Under  Diametral  Compression. 
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Table  VII 


Composition  of  Epoxy  Samples 


Sample 

Araldite  6020 
(pbw)1 

Araldite  508 

Araldite 

(pbw) 

(pbw) 

Epoxy  No.  3 

100 

100 

78 

Epoxy  No.  4 

100 

150 

50 

Epoxy  No.  5 

100 

50 

40 

Epoxy  No.  6 

100 

20 

Epoxy  No.  7 

100 

33 

33 

Parts  by  weight 


Table  VIII 

Composition  of  Polyester  Samples 


Sample 

Cadillac 

Laminae 

Laminae 

MEK 

EPX-126-3 

4134 

Peroxide 

(pbw) 

(pbw) 

(pbw) 

(drops/oz) 

Polyester  No.  8 

100 

3 

Polyester  No.  9 

100 

20 

7.5 

Polyester  No.  10 

100 

6 

Polyester  No.  11 

100 

30 

9 

Polyester  No.  12 

100 

30 

9 

Polyester  No.  13 

100 

20 

8 

Polyester  No.  14 

100 

20 

8 

Polyester  No.  15 

100 

10 

7 

Polyester  No.  16 

100 

10 

7 
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Casting  of  Flat  P ho toe las tic  Models 

J  number  of  flat  fiber  (or  rod)  reinforced  models  using  the  materials 
listed  were  cast.  Two  major  difficulties  encountered  in  the  castinq  process 
were  voids  (or  air  bubbles)  and  residual  stresses  in  the  model  9  P 

T+u  j°!ds  problem  was  solved  by  following  the  mixing  procedure  similar  tn 
the  method  devised  by  D'Agostino  et  al.  (Proc.  SESA  .m  2  955) 

The  procedure  involves  the  following  steps:  ’  ybb' 

After  weighing  (or  measuring)  the  resin  it  was  poured  into  a  rubber  sack- 
if  more  than  one  resin  was  used  they  were  mixed  together.  ’ 

TmmpJ?n/^ber  was,  suspended  until  the  bubbles  rose  out  of  the  resin, 
mm .rs mg  the  rubber  sack  in  hot  water  to  reduce  the  viscosity  of  the  resin  or 

Sf  S  davsnfnrVpnnym  t0  UP  the  Pr°CeSS  somewha^  It  took  a  minimum 

or  tour  days  for  epoxies  and  one  day  for  the  polyester  resin. 

syringed  harde"er  for  epoxy  was  smUarly  Prepared  in  advance  in  a  hypodermic 

The  rubber  sack  was  then  tied  below  the  liquid  level  and  hardpnpr  -in 
jected  into  the  resin  using  hypodermic  syringe  and  needle.  The  rubber  sack  was 

.mt?ia?hinhbe  °W  the  need1e  llole  and  mixing  was  done  by  kneading  the  rubber  sack 
until  the  homogeneous  mixture  was  obtained.  9  6  ruDDer  sack 

fnni“nS1'ierflelreduction  of  residual  stress  due  to  casting  was  obtained  bv 
following  Edelman  s  casting  procedures.  Edelman's  procedures  are: 

a.  Cool  the  assembled  mold  to  a  low  temperature  (approximately  30°F)  and 
store  at  that  temperature  for  24  hours. 

b.  Pour  the  resin  and  hardener  mixture  into  the  mold  and  cure  at  the 
low  temperature— for  24  hours. 

c.  Raise  the  temperature  at  a  rate  of  5°  per  hour  to  room  temperature. 

d’  tempera ture^aSt  ^  ^  mold  after  a  24"hour  cure  at  ™om 

A  Live-Load  Loading  Frame  Design 

is  Droceedinn9"  tL*!  J1*H®’1°5d  loadin9  frame  ha*  been  completed  and  fabrication 
proceeding.  The  loading-frame  as  designed  will  have  the  following  capabilities: 

a’  to ra up ^  nf^nrwT-6  /tensi°n  or  compression)  2,000  pound  and  maximum 
torque  of  3,000  inch-pound.  The  axial  force  and  torque  can  be  aDDliPd 

Actuator6ntirOT  dfemltane°USly  thr°Ugh  an  Axia1 -Force-Torque  Hydraulic 
b*  by X1 2 -i’nch^enjth?^1  WhiCh  Can  be  accon™dated  is  10-inch  diameter 
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a. 


c.  Model  can  be  rotated  about  two  axes  and  translated  alonq  two 
perpendicular  directions. 

d.  Two  directions  of  observation  are  possible. 

2.  Holography  Studies 

Holographic  techniques  provide  a  very  sensitive  means  of  measurinq  dis- 

3n^  stres^es-  Jn  the  field  of  Applied  Mechanics  this  opens  up  the 
possibility  of  experimental  means  to  directly  support  or  contradict  various 

i  Efforts  were  started  to  set  up  a  system  to  produce  holograms  which 
are  intended  to  supplement  the  scattered  light  effort. 

Optical  Fixtures:  To  hold  various  optical  elements  like  lenses, 
beam  splitters  etc.  optical  fixtures  were  designed  and  machined! 

Each  was  provided  with  various  motions  for  positioning  and  a  heavy 
base  for  stability.  y 

Holographic  Table:  This  is  a  heavy  base  on  which  the  entire  set- 
up  rests.  This  should  have  good  damping  characteristics  and 
rigidity  with  respect  to  bending  to  avoid  nodal  vibrations.  For 
these  reasons  a  steel  table  was  bypassed  in  favor  of  a  7'  x  4'  x  0.5' 
reinforcea  concrete  table  which  was  designed  and  fabricated. 

]ableSus£ension:  Various  means  of  supporting  the  1600  lb.  concrete 
table  to  isolate  it  from  building  vibrations  were  experimentally 
evaluated  These  included  heavy  steel  sprinas  and  cement  blocks 
followed  by  air  tubes.  Each  was  studied  for  frequency  response  and 
rinally  air  tube  suspension  was  found  adequate. 


c. 


d. 


Light  Source : 


_  For  reasons 
a  laser  is  the  most  suitabl 
laser  starting  from  the  50 
quality  optical  elements  is 
mounts  for  the  mirrors,  pri 
and  precision  alignment  of 
mounted  on  a  heavy  I  beam, 
laser  self-starting  feature 
in  progress,  with  the  help 
alignment  and  lasing  action 


of  intensity  and  purity  of  wavelength 
e  light  source.  The  assembly  of  a 
mw  plasma  tube  and  associated  high 
almost  complete.  Various  precision 
sms  and  different  clamps  for  support 
the  plasma  tube  have  been  machined  and 
A  six  kilovolt  power  supply  with  a 
has  also  been  assembled.  Work  is 
of  a  small  laser,  to  get  precision 


C. 


Dii.crete-Element_Microstress  Analysis  of  Unidirectional  Fiber  Comoosi tes 
(Professor  I.  p.  Kicher  and  Dr.  J.  F.  Stevenson] - - 


A  formulation  for  a  displacement-method  discrete  element  for  use  in  two- 
dimensional,  linear  microstress  analyses  of  unidirectional  fiber  comoosite 
materials  that  expedience  deformation  in  a  plane  normal  to  the  direction  of  the 
fibers  has  been  completed.  The  discrete  element  itself  is  capable  of  modeling 
a  rectangular  rt-gion  of  matrix  that  contains  the  entire  cross  section  of  a 
Tiber.  The  displacement  state  within  an  element  is  expressed  in  terms  of  nodal 
displacement  parameters  and  local  coordinate  functions  which  are  defined  tn 
satisfy  the  d  splacement-equi librium  equations  and  stress 'matching  conditions 
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at  the  fiber-matrix  interface.  Displacement  boundary  conditions  that  the  local 
coordinate  functions  must  satisfy  are  chosen  in  such  a  manner  that  inter-element 
displacement  continuity  at  the  common  edge  of  contiguous  elements  can  be  assured. 
Numerical  solutions  were  obtained  for  ten  problems  that  consist  of  assemblages 
of  one  or  more  of  the  following  plane  strain  elements:  square  element  with  a 
centered  circular  fiber,  square  element  with  an  off-cente.'ed  circular  fiber, 
square  element  with  no  fiber,  and  square  element  with  a  centered  circular  "hole" 
(i.e.,  a  highly  flexible  fiber).  The  solution  for  each  problem  is  generated  by 
direct  minimization  of  the  potential  energy  of  the  assemblage.  The  Fletcher- 
Reeves  algorithm  coupled  with  a  scaling  transformation  is  used  to  carry  out  the 
minimization. 

Details  of  this  work  can  be  found  in  tne  Ph.D.  thesis  of  J.  F.  Stevenson 
submitted  to  Case  Western  Reserve  University.  A  separate  topical  report  will 
not  be  issued  but  a  condensed  version  of  the  paper  has  been  submitted  for 
publication  in  an  archive  journal. 


D.  Application  of  ihe  Theory  of  Physically  Nonlinear  Elastic  Solids  to 
Composite  Materials 

(Professor  Sendeckyj ) 

The  investigation  of  stress  and  strain  concentration  factors  due  to  in¬ 
clusions  of  various  shapes  in  a  physically  nonlinear  matrix,  reported  in 
Reference  5,  has  been  continued  with  emphasis  being  placed  on  developing  an 
analysis  method  that  may  be  extended  to  nonlinear  solids  for  solving  linearly 
elastic  inclusion  problems.  Considerable  success  has  been  achieved  in  this 
respect.  The  general  solution  of  the  elastic  curvilinear  inclusion  problem 
has  been  found  in  the  case  of  antiplane,  or  longitudinal,  shear  deformation. 

This  result  has  betn  used  to  discuss  the  longitudinal  shear  behavior  of  fila¬ 
mentary  composites  in  References  5  and  16.  The  general  solution  has  also  been 
used  to  solve  the  s.itiplane  deformation  problem  of  a  filamentary  composite  with 
randomly  spaced  fibers  (summarized  in  Sections  1  and  2  below).  In  the  case  of 
plane  deformation,  the  general  functional  form  of  the  solution  for  an  elastic 
curvilinear  inclusion  has  been  found  (summarized  in  Section  3).  The  multiple 
circular  inclusion  problem  in  plane  elastostatics  is  presently  being  worked 
on,  with  limited  results  being  presented  in  the  next  chapter  (Section  V  E  of 
this  report).  Finally,  the  solution  methods  for  nonlinear  problems  have  been 
reconsidered  and  apparently  a  promising  approach  has  been  found  (Section  4 
below) . 

1 •  Multiple  Circular  Inclusion  Problems  in  Longitudinal  Shear  Deformation 

Consider  a  homogeneous  unbounded  matrix  containii,d  N  circular  cylindrical 
inclusions  with  confers  at  z  =  a.  (see  Figure  26).  Let  the  radius  of  the 
.th  .  1 

i —  inclusion  be  R^  and  assume  that  all  the  inclusions  have  the  same  shear 

modulus.  Let  subscripts  1  and  2  on  elastic  field  quantities  refer  to  the  matrix 
and  inclusions,  respectively.  Let  the  stress  state  at  infinity  be  a“z  =  1, 

ayz  =  It:  is  rec)uirecl  to  determine  the  elastic  fields  near  the  inclusion. 

Using  the  notation  of  Reference  5,  it  is  readily  established  that  the 
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complex  potentials  for  the  problem  are 


F,  (Z)  =  Z  -  t  K”  £m 
m=l  n=l 

C 

m,n 

(1) 

*-v.„ 

p2  ^  =  0+K)  {  z-  ?•  Km 

m=l 

Pm  C 

i  — !ML  } 

n-1  z"^m,n 

(2) 

where 

Pm  =  NfN-l)1"'1 

(3) 

C'.n  =  Rn  >  «l,n  =  a„ 

(  n  =  1,2 . N) 

(4) 

and,  in  general. 

cm,n  ”  ai  +  Ri  " 

M'2 

(5) 

Cm,n  Sn-l,k  Ri  ^m-l,k 

-  5,-r2 

(6) 

]f  jcm-l,k  ~ai  I  -Ri;  otherwise  the  terms  are  not  computed.  The  calculations 
indicated  by  (5)  and  (6)  are  to  be  carried  out  for  all  m;  k-1 ,a . P  ; 

n  and  i  -  1»2,...,N.  The  prime  on  the  summation  symbols  in  (2) 

Ref.  17.  The  stresses  may  be  computed  from  (1)  and  (2)  by  using  9 


zx 


la 


yz 


dF 

dz 


«a c*tx  sr 

the  positions  of  the  centers  of  the  inclusion];  while,  the  points  z 

(rn  >  2)  are  the  result  of  successive  aversions  of  z  =  e  n  with  respertto 

efficients1?"  Whl’CVS+the  operation  indicated  by  (5).  The  co- 

Cm,n  ^  -  2)  dre  Products  of  squares  of  distances  from  the  sinqular 

points  to  the  centers  of  the  inclusions  in  which  the  singularities  occur. 

2. IS.' tM.-S.-S 
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2.  Effective  Shear  Modulus 


In  a  fundamental  paper.  Hill  (18)  developed  exact  expressions  for  the 
effective  moduli  of  multiphase  composites.  In  the  case  under  consideration, 
the  pertinent  expression  may  be  written  as 


U  *  U-j  +  V  B(u2  “  U]  ) 


(7) 


where  y  is  the  longitudinal  shear  modulus;  v  is  the  concentration,  or  volume 
fraction,  of  the  fibers;  and  B  is  the  average  strain  concentration  in  the  fiber. 
The  use  of  (7)  requires  knowledge  of  the  average  strain  concentration  factor 
in  the  fiber,  which  may  be  obtained  from  (2)  as  follows:  the  stresses  within 
the  typical  inclusion  being  considered  are 


a  -la 
xz  yz 


(1+K)  {  1  +  z  K" 
n=l 


m 

E' 


m,n 


(8) 


Since  the  singularities  are  outside  the  inclusion,  the  stresses  can  be  expanded 
in  a  Laurent  series  about  the  center  of  the  inclusion  to  give 


m  C  9 

°xz  -  %2  3  (1+K)  {  1  +  ”  Km  I'  [1  +  |S-  +  . 

y  m=l  n=l  cZ  cm,n 

m,n 

The  terms  with  2  yield  stresses  that  depend  on  cosine  and  sine  of  ke.  Upon 
averaging,  these  give  a  zero  contribution.  Thus,  we  get 


'xz  -  1  v  =  (1+k)  (  '  +  K"  r; 

m= I  n=  I 


m 


(9) 


for  the  averaged  complex  stress  combination. 

For  a  regular  array  of  fibers,  it  can  be  shown  after  some  rather  tedious 
computations  that 


xz 


(1+K)  { 


KS, 


3K2S^ 


(1-KS2)2  [(4a2)2  -  S2  K2] 


(cont'd) 
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+ 


(1-KS2)Z  [(4a‘i)4  -  4  K2] 


7K2S2 

(l-KS2)^(4a2)6  -  S2  K2] 

where 

Sn  =  an/4a2 


=  (1+K)  A  (10) 


(11) 


and 


a2k 


f1  1 

n=-»  [rrH-ncelaJ 

m=-« 


(12) 


where  c  and  a  are  parameters  defining  the  regular  array.  These  are  defined  in 
Figure  27.  The  prime  denotes  that  the  summation  is  to  be  carried  out  for  all 
integer  values  of  n  and  m  with  the  exclusion  of  n=m=0. 


Using  (1)  to  estimate  the  average  strain  concentration  factor  for  the 
fibers,  we  get 

u  =  U1  (1  +  2v  KA)  (13) 

for  the  effective  shear  modulus.  This  expression  yields  results  in  excellent 
agreement  with  numerically  obtained  graphs  given  in  References  19  and  20. 

3.  Plane  Deformation  -  Curvilinear  Inclusion  Problem 

Some  preliminary  results  on  the  curvilinear  inclusion  problem  in  plane 
elastostatics  were  presented  in  Reference  5.  The  present  section  deals  briefly 
with  results  generated  since  Reference  5.  Since  the  notation  that  is  used  can 
be  found  in  any  modern  text  on  the  theory  of  elasticity,  it  will  not  be  given 
here. 

Consider  an  elastic  matrix  containing  a  single  curvilinear  inclusion  whose 
boundary  is  given  parametrically  by 

z  =  w(c)  »  |t|  =  1.  (14) 

Letting  subscripts  1  and  2  on  elastic  quantities  refer  to  the  matrix  and  in¬ 
clusion  respectively,  the  boundary  conditions  corresponding  to  a  perfect  bond 
at  the  interface  are 

U1  +  ivi  =  u2  +  iv2 

for  |c|  =  1  .  (15) 

°ppl  1ap9l  app2  1ape2 
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Figure  27.  Regular  Array  of  Circular  Fibers 
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The  complex  potentials  <f>(c)  and  *(c)  satisfying  (15)  are 


<Mz)  -  ♦!  (C)  -  f(z)  -  ^  ?  [0(c)] 

-  frf R  *  U6) 

¥i  (z>  a  *1  <«>  -  f k  f  [ajU)] 

+  f|  IHcT  3c  fi4cf  ?  [n(?)] 

+  f^f  f  +  f?f  R  Cfi(c)3  ,  (17) 

*2^  s  Tip  *  (18) 


Vz)  =  T+3  h(z^  *  (19) 

where 

^  =  Cl  (kt+1  }  -  (*<2+l)]/[r(k1+l)  +  k2+1]  , 

3  *  Cr (<-| ~l )  -  (<2-l  )3/Lr(<^+l )  +  <2+l]  »  (20) 

r =  ^^i 

and 

fi(?)  =  oj  ( 1/ c )  .  (21) 

General  solution  (16-19)  can  be  readily  specialized  for  particular  com¬ 
binations  of  material  properties.  For  a  cavity  or  an  inclusion  of  vanishing 
rigidity,  one  gets 

4>U)  a  f (c)  -  R(l/c)  , 

*<c)  ■  h({)  -  t  Ij-  BO/O  -  run  •  (22) 

For  a  rigid  curvilinear  inclusion,  r  ■  »  and  the  general  solution  reduces 

tu 
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♦  (?)  •  f(c)  *  7  fi  (1/?)  > 

▼(e)  ■  h(e)  +  f(i/{)  -^fj'Crajl'M/t)  +  f'(t)  ]  .  (23) 

A  more  detailed  treatment  of  these  results  can  be  found  in  Reference  21. 

Some  preliminary  results  on  the  application  of  the  general  solution  to 
composites  are  reported  in  the  section  by  G.  P.  Sendeckyj  and  I.  W.  Yu  in  this 
report. 

4.  Nonlinear  Antiplane  Deformation 

During  the  period  covered  by  this  report,  four  methods  of  solving  physic¬ 
ally  nonlinear  elasticity  problems  were  examined,  namely, 

1)  successive  approximations, 

2)  perturbation  series, 

3)  series  in  terms  of  trigonometric  functions,  and 

4)  series  in  terms  of  inverse  powers  of  the  radial  coordinate.  Of 
these  methods,  the  fourth  seems  most  promising;  but  only  partial  results  are 
available  at  this  time.  Of  the  other  three  methods,  expansion  of  the  dis¬ 
placements  in  terms  of  trigonometric  series  in  the  0  direction,  with  unknown 
coefficients  depending  on  the  radial  variable,  leads  to  an  infinite  system 

of  badly  coupled  nor.iinear  differential  equations.which  cannot  be  solved.  The 
successive  approximations  technique  is  extremely  tedious  to  use,  while  the 
perturbation  series  approach  is  not  much  batter.  Future  work  in  this  area 
will  be  centered  on  use  of  the  fourth  method. 

E.  Multiple  Circular  Inclusion  Problem  in  Plane  Elastostatics 

(Professor  G.  P.  Sendeckyj  and  Mr.  I.  W.  Yu) 

Recent  technical  literature  abounds  with  solutions  of  problems  of  plane 
elastostatics  for  regions  with  circular  boundaries.  Yet,  upon  scanning  the 
literature,  it  is  readily  seen  that  very  few  solutions  are  available  for  the 
case  of  more  than  one  elastic  inclusion.  These  solutions  are  restricted  to 
special  geometric  configurations  of  the  inclusions.  The  present  work  is  an 
attempt  to  generate  analytical  solutions  for  the  case  of  randomly  spaced 
elastic  circular  inclusion  in  a  solid  undergoing  plane  deformation.  The 
approach  used  to  solve  the  corresponding  problem  for  antiplane  deformation  (see 
Section  V  D  in  this  report)  is  being  used.  At  present  only  limited  preliminary 
results  are  available.  These  are  summarized  below. 

1 .  Two  Circular  Cavities  on  x-axis 

Consider  an  elastic  slice  weakened  by  two  circular  cavities  of  unit  radius. 
Let  the  centers  of  the  cavities  be  put  at  the  points  xs  +  a,ys0(a>  1). 

Let  the  applied  stresses  at  infinity  be  specified  as 


-  66  - 


®  «  00  .00 

°xx  =  -1’  °yy  =  *  axy  =  °* 


(1) 


It  is  required  to  determine  the  elastic  fields  in  the  vicinity  of  the  cavities. 

correS^„li„TtotO)aar?en0tat1'°l!  °?  e,aStiCity  the°ry’  the  c°mp,ex 


K  <z)  =  0.  ^(z)  =  z  . 


(2) 


Upon  using  (2)  as  the  input  to  the  general  form  of  the  solution  of  a  circular 
cavity  (see  Section  V  D  in  this  report),  it  follows  that  circular 


♦  (z)  =  -  -L  +  _L 

z-a  z+a  5 

Hz)  =  z  -  a  [— 


(3) 


C^+^7]-  (4) 

Equations  (3)  and  (4)  constitute  the  first  approximation  of  the  desired  solu- 
t  on  The  stresses  due  to  these  potentials  satisfy  the  boundary  conditions 

anriy(aiPh0Xlm?te  rh  A"  ""p™ved  approximation  follows  upon  connecting  (3) 
can  be’writuHs  9e"era  f°rm  °f  the  S°,Ution-  The  PK°"d  aPProxi™Jion 


*(2)  =  ■  w  '  w  f  rX  •  4-  3  c  +  TrSFTO  ] 


+  a  3 

+  "X  +  3" 


jo  bo 
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b0 


be  b0  (z-a+l/bo)2  (z+a+l/b0)2 
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- *-  +  - : - - 
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Hz)  =  z  -  a  — L-* - 1  .  _1  , 
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~2~  - - - 
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(cont'd) 


-  67  - 


1 
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1 


(z+a-l/boV 


(z-a+l/b0)‘ 


1 

(z-a+l/bo)3 
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1 

(z+a-l/bo)3 


+  _3  __J _  _  1 

2bo  (z-a+l/bo)4  (z+a-l/bo)4 

where 

b0  =  2a  . 


(6) 


(7) 


The  third  approximation  was  also  generated;  but  the  structure  of  the 
general  solution  is  still  not  clear.  A  fourth  approximation  is  being  gene¬ 
rated  at  the  present  time.  This  should  make  it  possible  to  guess  the  struc¬ 
ture  o>  the  general  solution.  Once  this  solution  is  known,  generalization 
to  a  randomly  perforated  slice  is  expected  to  be  straightforward. 

F.  Structural  Analysis  for  Laminated  Fiber  Composite  Structures 
(Prof.  L.  A.  Schmit  and  Mr.  S.  N.  Patnaik) 

Laminated  fiber  composite  plates  and  shells  composed  of  a  collection  of 
highly  orthotropic  lamina  with  noncoincident  laminate  axis  make  the  total 
structure  anisotropic  and  heterogeneous  through  the  thickness.  Current  analy¬ 
sis  of  fiber  composite  plates  and  shells  are  based  on  several  basic  simplifying 
assumptions  taken  from  isotropic  plate  and  shell  theory.  Since  the  behavior 
of  isotropic  and  heterogeneous  anisotropic  materials  differ  to  a  great  extent, 
it  is  but  logical  to  question  the  validity  of  the  isotropic  plate  and  shell 
assumptions  where  applied  to  laminated  systems. 


The  goal  of  the  present  research  is  to  assess  the  validity  of  the  basic 
assumptions  presently  used  in  the  analysis  of  laminated  plates  and  shells  and 
to  propose  an  improved  set  of  hypothesis.  The  structure  is  discretized  by 
three-dimensional  elements.  An  element  stiffness  matrix  with  linear  aniso¬ 
tropic  material  properties  is  generated  using  the  theorems  of  potential  energy 
and  trivariant  linear  Hermite  displacement  functions  developed.  Since  the 
thickness  of  laminated  structures  is  usually  small  compared  to  the  other  two 
dimensions,  to  avoid  ill  conditioning  of  the  stiffness  matrix  the  aspect 
ratio  has  to  be  kept  small.  That  is,  the  mesh  points  of  the  finite  element 
grid  must  be  close  together  and  hence  linear  interpolation  for  displacements 
would  serve  the  purpose.  This  results  in  24  degrees  of  freedom  per  element. 


Basic  stiffness  matrix  is  generated  and  numerical  verification  is  com¬ 
pleted  using  various  simple  test  cases.  To  verify  the  extent  to  which  strain 
variations  th  ough  the  thickness  of  laminated  structures  deviate  from  the 
usual  linear  assumption,  a  beam  shown  in  Figure  28  with  four  layers  through 
the  thickness  is  analyzed. 


Figure  29  shows  that  along  x-direction  the  displacement  is  linear,  but 
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Figure  29.  In  Plane  Displacement  Through  Thickness 


the  displacements  along  y-direction  are  nonlinear.  The  displacement  in 
y-direction  suggest  that  there  is  an  inplane  unbalanced  shear  force,  and 
this  staggers  the  laminas.  This  is  due  to  the  fact  that  the  stiffness  alonq 
the  y-direction  is  much  less  and  antielastic  curvature  cannot  be  represented 
adequately  by  the  model,  as  we  have  only  one  element  along  the  v-direction. 
The  analysis  of  a  plate  is  underway. 


G.  Three-Dimensional  Finite  Element  Analysis  -  Mixed  Stress-DisDlacement 

Approach - - c - 

(Professor  Schmit  and  Mr.  P.  Obaid) 


This  report  contains  the  analytical  investigation  carried  out  in  order 
to  examine  the  validity  of  basic  assumptions  upon  which  the  current  analysis 
methods  for  laminated  fiber  composite  plates  and  shells  rest.  Laminated  fiber 
composite  plate  and  shell  configurations  may  be  viewed  as  a  collection  of 
highly  orthotropic  lamina  and,  in  general,  assuming  the  configuration  is  not 
unidirectional,  the  structure  is  highly  heterogeneous  through  the  thickness. 
The  stress-strain  characteristics  of  the  orthotropic  lamina' tend  to  be  linear 
in  the  fiber  direction,  somewhat  nonlinear  transverse  to  the  fibers,  and 
highly  nonlinear  under  inplane  shear.  A  usefuladvance  in  behavior  predic¬ 
tions  might  be  achieved  by  introducing  shear  nonlinearity  while  retaining  the 
assumption  of  linear  stress-strain  behavior  in  the  longitudinal  and  trans¬ 
verse  directions. 

After  assessing  the  validity  of  basic  assumptions  of  laminated  plate  and 
shell  analyses,  the  task  will  be  set  forth  to  find  an  improved  hypothesis  and 
to  generate  the  general  capabilities  so  as  to  analyze  such  plate' and  shell 
configurations. 

Analysis  is  based  on  the  mixed  stress-displacement  approach.  This  formu¬ 
lation  corresponds  to  the  Reissner  energy  principal  (Ref.  22).  Reissner's 
variational  principal  states  that  of  all  possible  states  of  stress,  that  which 
makes  the  Reissner  energy  stationary  satisfies  the  equilibrium  equations  and 
force  displacement  relations  as  well  as  the  natural  boundary  conditions  (both 
displacement  and  force)  and  is  therefore  the  actual  state  of  stress  and  dis¬ 
placement. 


Formulation  of  a  three-dimensional  rectangular  parallelepiped  element 
is  pursued  to  begin  with.  The  Reissner  energy  jid  of  a  typical  interior 
discrete  element  is  given  as  K 
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where  assumed  nodal  displacements  and  stresses  consist  of  triple  products  of 
zeroth  order  Hermite  interpolation  polynomials. 


u(x,y,z)  -  E  E  E  H^'(x)  H^'(y)  H^'(z)  u... 

i*i  j=l  k=l  01  0J  ok  ^ 

a  (x,y,z)  =  E  E  E  H<°fx)  H^y)  H^}(z)  a 

X  i-1  j=l  k-1  01  0J  ok  xijk 


u ( °)t„\  u(o), 


(2) 


(3) 


In  the  above  equations  u  and  a  represent  the  displacement  and  stress  res- 

X  XL 

pectively  in  the  x  direction.  ZeroLn  order  Hermitian  polynomials  are  de¬ 
fined  as 


H„,  (x)  -  1  -  f 


Ho2  (x) 


f-  etc. 

a 


(4) 

(5) 


a  being  the  distance  between  two  interpolation  stations,  u...  and  a 

1Jk  xijk 

represent  the  nodal  displacement  and  nodal  stress  along  x-axis. 

For  a  general  homogeneous  anisotropic  elastic  body,  stress-strain  rela¬ 
tions  including  thermal  contributions  can  be  written  as 
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where  t  and  a  represent  the  strain  and  stress  vectors  respectively,  a 
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contains  the  thermal  contribution  and 
coefficients  out  of  which  only  21  are 
[b]. 


[b]  contains  the  generalized  36  elastic 
independent  due  to  the  symmetry  of  matrix 


f coefficients  or  so-called  stiffness  coefficients  are  reduced 
fr°m  *■  b0.9  inciePenc,ent  constants  for  an  orthotropic  body  and  to  only  2 
independent  constants  for  the  case  of  an  isotropic  body.  These  elastic  con¬ 
stants  are  given  in  reference  (23).  However,  it  must  be  recognized  that 
3-d  sh?anng  stresses  are  expressed  in  different  order  in 
equation  (6)  than  in  reference  23)  due  to  which  position  of  various  rows  and 
columns  In  matrix  [b]  and  in  vector  J  are  also  changed.  Hencl  to  awid  » 
confusion  the  expressions  for  various  elastic  coefficients  are  given  as  V 
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where  Exx>  Eyy  and  Ezz  are  Young's  moduli  corresponding  to  X,  Y  and  Z  axes 

respectively.  Gxyn  G^z  and  Gzx  are  shear  moduli  corresponding  to  planes  xy, 

yz  and  zx  respectively.  vyx,  vxy3  .....  vzx  are  the  Poisson's  coefficients 

which  characterize  transverse  compression  or  tension  in  the  direction  of 
coordinate  axes.  Coefficients  of  mutual  influence  n  .  n  n 

X » xy  *  *  *  *  *  *  2^ 

are  shears  in  the  planes  parallel  to  the  coordinates  under  the  influence  of 
normal  stresses.  Chentsov's  coefficients  p  ,  p  ...  u  charac 

yz,xy  ^xy.yz . zx  cnarac 

terize  the  shear  in  the  plane  which  are  parallel  to  the  coordinates’and  which 
include  the  tangential  stresses  parallel  to  other  coordinate  planes. 

Complementary  strain  energy  density  uc  is  given  as 
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Substituting  equation  (6)  into  equation  (9)  we  get 
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By  substituting  equations  (2),  (3)  and  (10)  into  equation  (1),  the  fol¬ 
lowing  expression  is  obtained  after  integration  over  the  volume  of  a  discrete 
element. 


where 


nR  = 

t  [CJ  a  - 

1  -T 

J  a 

[D] 

-*■ 

0 

T 

+  0 

1x24  24x48  48x1 

c  1x48 

48x48 

48x1 

1x43 

stT 

=  L  Uin,  U]12,  U121, 

U122’  U21 1 

*  U21 2  * 

U221  * 

U222  *  v 

. v222 *  wlll 

’  W1 1 2  *  * 

•  •  • 

’  w222 

3  ...  (11) 


v- 


112’  '  *  * 
(12) 


-  74  - 


Nil'  °x  1 T  2  *  "*121 '  °x122’  %„■  %12’  °x221  *  °x222'  "y, , ,  * 


v,,  •  *  *  9  9  i  •  •  •  *  i  )  cr  . 

112  y222  2m  ZH2  2ooo* 


Txylll  ’  Txyl  1 2 . Txy 222 *  Tyzin’  T^zh2 . 

yz222 *  Tzxm’  Tzx112 . Tzx222  ^  (13) 

and  define  column  vectors  P  (x)  ,  *'(x).  fry),  $'(yh  ft(z)  and  ft'(z)  as 

fT(x)  ■  [H<?)  ,x),  H<°>(x),  H<?)(x).  H<?)(x;,  h!»)(x),  H<°>(x),  h<°>(x}h<<>>(x,] 
5T(y)  ■  H<°>(y),  $>(y).  H<°>(y),  H<°>(y),  H<°>(y),  „(°)(y),  H<°|y)] 

rT<2)  ■  CH<^(Z),  h<°>(2),  h<°>(2),  hJ®)(*).  h<°>(2),  h<°>(2),  h<°>(2),  H<»)(2)] 

•  •  •  (14) 

P'(x).  f  ' (y)  and  ft1  (z)  are  the  first  derivatives  of  P(x)  fi(y)  an 
respectively.  Now  displacements  and  stresses  can  be  expressed  in  the  form 
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where  subscripts  i,  i-8  etc.  represent  the  particular  elements  of  vectors 
Q  and  R.  Matrix  [c]  is  partitioned  as  follows 
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where  all  of  the  submatrices  are  8x8  dimensional  arrays.  Elements  of  each 
submatrix  are  given  as 
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The  range  of  subscripts  i  and  j  is  from  1  to  8.  Now  [D]  can  be  partitioned  as 
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[D]  is  a  symmetric  matrix  and  all  the  submatrices  such  as  et 

are  8x8  dimensional  arrays.  The  elements  of  these  matrices  can  be  defined 
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where  subscripts  i  and  j  both  vary  from  1  to  8.  In  equation  (19)  b^  represents 
the  stiffness  coefficients  which  have  already  been  defined  in  equation  (8). 


The  column  vector  Q  is  the  thermal  contribution  and  is  given  in  the  form 
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(21) 


Range  of  subscripts  is  from  1  to  8. 

A  simple  three-dimensional  elasticity  problem  was  solved  using  the  Reissner 
variational  principal.  The  test  problem  consisted  of  a  prismatic  bar  stretched 
by  an  axial  force  P  applied  at  one  end.  Dimensions  of  the  bar  and  other  neces¬ 
sary  details  are  given  in  Figure  30.  One  end  of  the  bar  was  fixed.  Reissner 
energy  for  this  problem  is  given  as 

nR  «  t  [C]  t  -  1  f  [D]  o  +  ST  0  -  xT  i*  (22) 

where  P  is  a  vector  containing  24  generalized  loads  associated  with  each  of  che 
24  nodal  degress  of  freedom  contained  in  X.  Other  notations  have  been  defined 
earlier. 

Take  the  variation  of  Reissner  energy  nR  (Refs.  24  and  25)  and  define: 

(11)  (12)  (21)  T  (22) 

[R]  ■  -  [D]  ,  [R]  -  [C]  ,  [R]  =  [C]T  and  [R]  =  [0]  ..  (23) 

Hence 
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R(ll)^  r(12)^  r^I)^  and  r(22)  are  caiie(j  as  Reissner  submatrices.  Equation 
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(24)  was  successfully  solved  using  computer  technique  and  results  showed  ex¬ 
cellent  agreement  with  the  theoretical  analysis  (Ref. 23)  for  isotropic,  ortho¬ 
tropic  and  anisotropic  cases. 

A  prismatic  bar  was  then  divided  into  four  elements  in  the  xy  plane  as 
indicated  in  Figure  31.  Element  Reissner  matrices  were  combined  with  the  help 
of  a  general  variable  correlation  table  (Refs.  27  and  28)  to  give  the  master 
Reissner  matrix  for  the  system.  Boundary  conditions  for  this  case  were  sliqhtly 
altered.  Neither  end  was  kept  fixed  and  instead  of  prescribing  the  load  at  one 
end,  uniform  half  unit  displacements  were  prescribed  at  both  ends  in  the  oppo¬ 
site  directions.  Horizontal  constraints  were  applied  for  displacements  u  and 
v  at  one  end  such  that  along  the  plane  z  =  0,  the  displacement  u  is  zero  if 
x  is  zero  and  displacement  v  is  zero  if  y  is  zero.  The  final  analysis  qives 
the  rest  of  the  u  and  v  displacements  and  all  the  stresses  and  applied  loads 
at  each  node.  Nodal  loads  thus  obtained  come  out  to  be  equal  to  the  work  equi¬ 
valent  loads.  The  solution  was  checked  with  the  solution  given  in  reference 
23  and  complete  agreement  exists  between  the  finite  element  analysis  and  the 
theoretical  analysis  for  isotropic,  orthotropic  and  anisotropic  cases. 

As  the  number  of  elements  increases,  the  total  number  of  unknown  vari¬ 
ables  also  increases.  For  the  case  of  one  element,  the  Reissner  matrix  was 
of  the  size  72  x  72  and  for  four  elements  its  size  increases  to  162  x  162. 
Hence,  if  the  number  of  elements  increases  further, the  size  of  the  master 
Reissner  matrix  becomes  very  large  ana  the  problem  arises  as  to  how  to  deal 
with  such  a  large  matrix  in  an  economical  way.  It  is  also  desirable  to  make 
the  computer  program  fairly  general  and  independent  of  the  total  number  of 
elements.  In  these  connections  the  following  points  need  to  be  considered: 

(1)  The  Master  Reissner  matrix  is  symmetric;  hence,  only  the  upper 
triangular  elements  need  to  be  stored.  This  will  cut  down  the 
storage  requirement  to  a  great  extent. 

(2)  The  Master  Reissner  matrix  has  a  7arge  number  of  elements  which 
are  zero.  Hence, an  efficient  subroutine  should  be  used  for 
solving  the  large  number  of  simulcaneous  equations.  In  this 
regard,  it  will  be  worthwhile  to  make  use  of  the  wave  front 
technique  discussed  in  reference  (28). 

(3)  Maximum  possible  use  will  be  made  of  extra  storage  capacity  of 
the  computer  in  order  to  handle  very  large  matrices. 

Further  effort  is  now  directed  toward  writing  a  general  program  based  on 
the  above  considerations,  which  will  facilitate  assessment  of  the  basis 
assumptions  on  which  current  laminated  plate  and  shell  theory  rests. 
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Figure  31.  Prismatic  Bar  Divided  into  Four  Finite  fclements. 
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SECTION  VI 


FAILURE  OF  COMPOSITE  STRUCTURAL  ELEMENTS 

A.  Interfiber  Failure  of  Composite  Materials 

(Professor  T.  P.  Richer  and"  Dr.  J.  8."  Koeneman) 

The  conditions  which  exist  at  failure  of  unidirectional,  continuous  fiber, 
off-axis  epoxy  matrix  composite  materials  were  investigated.  A  micromechanic 
failure  criterion  was  established  which  reflected  these  conditions,  and  the 
results  were  used  to  explain  and  predict  composite  behavior.  Observations 
through  a  microscope  were  made  of  interface  failure  of  single  fiber  inclusion 
specimens  during  loading.  The  analysis  of  creep  tests  was  used  to  determine 
the  point  at  which  the  matrix  went  into  the  bulk  plastic  region,  and  this  was 
used  to  give  an  indication  of  the  stress  state  in  the  matrix.  Based  on  these 
results  and  a  comparison  of  polymer  failure  characteristics,  the  failure  criterion 
chosen  was  a  critical  creep  rate  normal  to  the  plane  of  fracture.  The  Eyrinq 
rate  process  theory  in  conjunction  with  a  combined  stress,  nonlinear  creep 
theory  was  used  to  apply  the  theory  to  composite  materials. 

To  calculate  the  plastic  stress  state,  a  procedure  proposed  by  Neuber  was 
used  which  relates  the  equivalent  plastic  stress  to  the  elastic  solution  by 
means  of  a  leading  function  which  is  determined  by  the  nonlinear  stress-strain 
law.  This  method  was  verified  on  several  test  problems.  Applications  of  the 
failure  criterion  indicates  that  the  transverse  strength  of  unidirectional  com¬ 
posites  should  increase  as  the  strain  hardening  of  this  type  matrix  increases. 

The  effect  of  residual  stresses  on  failure  strength  is  predicted  to  be  minimal. 

The  shear  stresses  existing  around  fiber  ends  or  where  a  fiber  passes  through 
a  void  have  a  significant  effect  on  matrix  yield  and  composite  failure.  The 
difference  in  temperature  dependence  of  crosspiied  strength  between  glass- 
fiber  and  graphite-fiber  composites  is  explained.  The  width  effect  of  cross- 
plied  strength  is  explained,  and  a  width  effect  and  layer  effect  on  crossplied 
creep  rate  was  predicted  and  then  experimentally  verified.  An  explanation  is 
provided  for  the  strengthening  effect  on  epoxies  of  elastomer  inclusions.  The 
sensitivity  of  graphite-fiber  strength  to  the  stress  state  in  the  matrix  is 
predicted.  A  method  of  determining  the  yield  stress  of  high  strain  hardening 
epoxies  by  the  analysis  of  creep  data  was  developed.  This  technique  was  also 
adapted  to  measure  the  bulk  yield  of  the  matrix  in  a  unidirectional  composite. 

This  work  has  been  presented  in  detail  in  a  Topical  Report  (29).  Copies 
are  available  from  the  Department  of  Solid  Mechanics,  Structures  and  Mechanical 
Design,  Case  Western  Reserve  University. 


B.  Fracture  Toughness  of  Fiber  Composites 

(Professor  T.  P.  Richer  and  Mr.  T.  C.  Esselman) 

In  the  past  year,  many  tests  have  been  conducted  to  determine  the  nature 
of  the  stress  intensity  factor  and  the  plane  strain  fracture  toughness  of  fiber 
composites.  An  effort  has  Deen  made  to  adapt  the  ASTM  Standard  (30)  for 
fracture  toughness  testing  of  metals  to  an  orthotropic  or  anisotropic  body. 

The  specimens  are  tested  in  three  point  bending  (Figure  32).  A  chevron  notch 
is  machined  into  the  specimens  and  a  crack  is  fatigued  in  this  notch.  A 
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solution  of  black  India  in»  and  ethylene  glycol  is  injected  into  the  crack  as 
it  is  being  fatigued  so  that  its  progress  may  be  followed  visually  and  so  that, 
after  fracture,  the  initial  crack  length  can  be  easily  measured  by  the  exten¬ 
sion  of  the  ink  solution  (Figure  33).  Since  the  ink  is  allowed  to  dry  before 
the  test,  this  assumes  that  there  is  no  slow  crack  growth  during  the  prefracture 
loading.  Visual  observations  during  the  test  justify  this  assumption. 

During  the  test,  a  continuous  record  is  kept  on  an  x-y  recorder  of  the 
applied  load  and  the  separation  displacement,  fhe  displacement  apart  at  the 
sides  of  the  crack.  The  separation  disclacement  is  measured  using  a  double 
cantilever  beam  transducer  with  strain  gages.  Knife  edges  are  clamped  to 
both  sides  of  the  crack  and  the  transducer  is  attached  to  these.  The  load  is 
measured  with  a  load  cell.  These  transducers  can  be  seen  in  position  in 
Figure  32. 

The  tests  were  conducted  on  an  Instron  at  a  loading  rate  of  .100  inches 
per  minute.  The  strain  rate  was  varied  only  for  a  series  of  tests  where 
effect  of  variation  with  strain  rate  was  of  interest.  The  Instron's  compres¬ 
sion  load  cell  was  also  used  as  a  check  of  the  load  measured  on  the  x-y 
recorder. 


Both  unidirectional  and  crossplied  fiber  glass  were  tested. 


1.  Unidirectional  specimens 


The  unidirectional  specimens  were  chosen  with  the  crack  running  parallel 
to  the  fibers  since  it  has  been  shown  by  a  number  of  authors  (31,32)  that  the 
value  of  the  stress  intensity  factor  for  an  orthotropic  plate  with  the  crack 
propagating  in  one  of  the  material  principal  directions  can  be  found  by  using 
the  isotropic  analysis.  The  isotropic  analysis  used  by  the  ASTM  is,  for  the 
stress  intensity  factor, 


K 


I 


3  P  L  a^ 
B  W 2 


for  the  three  point  bend  test.  Y  is  found  by  boundary  collocation  (33)  and 
is  given  by 

y  =  1.93  -  3.07(a/W)  +  14.53(a/W)2  -  25.11(a/W)3  -  25.80(a/W)4 

In  order  to  make  the  tests  completely  meaningful,  a  state  of  plane  strain 
should  exist  in  the  specimen.  However,  several  things  indicate  that  for  speci¬ 
men  thicknesses  less  than  one  inch  (the  maximum  thickness  tested  in  these 
experiments)  a  state  of  plane  strain  is  not  present.  The  load-displacement 
curves  show  this.  Figure  34  is  an  example  of  a  load-displacement  curve.  This 
curve  was  chosen  as  an  example  since  the  change  in  slope  was  so  pronounced. 

In  most  others,  the  change  in  slope  is  more  subtle  but  nevertheless,  present. 

The  nominal  stress  at  the  crack  tip  was  found  using  a  simple  relation 
which  gives  1 

„  =  3  P  L 

N  B  (W-a)2 
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Figure  33.  Photo  of  Fractured  Specimen. 
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For  the  specimen  whose  load-displacement  curve  is  given  in  Figure  34,  the 
nominal  stress  at  fracture  was  6154  psi  and  the  nominal  stress  at  the  break 
in  slope  is  approximately  1170  psi.  It  has  been  found  by  Koenemari  (29)  that 
the  yield  stress  for  unidirectional  fiber  glass  composites  with  the  fibers  per¬ 
pendicular  to  the  tensile  loading  direction  is  1800  psi.  Considering  that  the 
nominal  stress  Includes  no  stress  concentration  factors,  the  break  In  the  slooe 
corresponds  to  yielding  in  the  material.  The  region  above  the  yield  point, 
which  is  fairly  linear  over  most  of  Its  range,  can  be  attributed  to  either  slow 
crack  growth  or  gross  plastic  deformation.  If  it  is  gross  plastic  deformation 
the  tip  of  the  crack  would  be  blunted  by  the  plastic  zone. 

A  plot  of  the  field  stress  vs.  the  crack  depth,  shown  in  Figure  35,  also 
Indicates  that  plane  stress  does  not  exist  In  the  specimen.  ,  The  basis  of  the 
linear  elastic  plane  strain  fracture  toughness  Is  that  of(a)‘  Is  constant,  where 

Of  Is  the  stress  In  the  field  without  the  crack  and  a  Is  the  crack  length.  Thus, 

a  plot  of  log  of  vs.  log  a  should  have  a  slope  of  -.5.  Figure  35  shows  a  slope 

of  -1.44  which  Indicates  that  linear  elasticity  does  not  apply.  The  nominal 
stress  at  the  crack  tip  is  plotted  in  Figure  36  as  a  function  of  crack  depth, 
and  the  value,  of  K  as  a  function  of  crack  depth  is  plotted  In  Figure  37. 

The  value  of  K  vs.  the  thickness  of  the  specimen  is  shown  In  Figure  38. 

Each  point  on  the  graph  represents  a  series  of  six  to  fifteen  specimens  which 
were  normalized  with  respect  to  the  crack  depth  ratio.  Therefore,  the  value 
of  K  at  each  thickness  corresponds  to  the  value  of  K  at  a  constant  crack 
depth  ratio,  chosen  arbitrarily  as  5.  The  trend  that  this  curve  takes  Is 
opposite  to  what  is  expected  from  metals.  The  plane  strain  fracture  toughness 
is  generally  accepted  as  a  conservative  value,  l.e.,  it  is  lower  than  the 
stress  intensity  at  a  state  of  plane  stress.  Therefore,  the  slope  in  the  tran¬ 
sition  region  from  plane  stress  (thin  specimens)  to  plane  strain  (thick  speci¬ 
mens)  should  be  negative.  The  reason  for  this  discrepancy  has  not  been  found 
yet  but  Is  probably  in  the  nature  of  the  plastic  zones  for  composite  specimens. 

The  variation  of  K  with  respect  to  the  rate  of  loading  was  also  studied. 
The  results,  plotted  on  semi-log  paper,  are  shown  in  Figure  39.  The  trend  is 
as  it  was  expected  to  appear. 

The  variation  of  K  with  respect  to  void  content  of  the  material  is  also 
being  studied  but  only  preliminary  results  are  now  available.  This  data  will 
appear  in  a  later  report. 

The  data  presented  for  the  unidirectional  specimens  should  not  be  taken 
as  values  of  plane  strain  fracture  toughness.  As  was  pointed  out,  the  speci- 
mnswere  in  a  transition  region  between  plane  stress  and  plane  strain.  The 
ASTM  suggests  for  metals  that  specimen  thickness  should  be  greater  than 

2.5(Kj/  o^)2  for  a  condition  of  plane  strain  to  exist.  Assuming  this  to  be 

true  for  composites  also  would  indicate  that  specimen  thickness  of  approxi¬ 
mately  2.5  inches  or  greater  are  necessary.  The  difficulties  in  fabricating 
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Figure  34.  Load -Displacement  Record  for  Specimen  0-1210-E. 
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and  testing  a  specimen  this  large  make  it  prohibitive.  Even  though  the  data 
presented  is  not  immediately  useful  in  design,  it  may  be  helpful  in  under¬ 
standing  the  phenomena  of  the  plane  stress  fracture. 

2.  Crossplied  specimens 

Preliminary  investigations  have  been  done  on  crossplied  specimens  of 
varying  fiber  orientation.  Figure  40  shows  the  plot  of  K  vs.  fiber  angle 
orientation  a  is  the  plus  and  minus  fiber  angle  orientations  from  the  di¬ 
rection  of  crack  propagation  (see  Figure  40  insert).  Isotropic  stress 
intensity  factor  analysis  has  beer  used  although  theory  shows  this  not  to  be 
valid  for  the  general  anisotropic  case.  Further  studies  similar  to  those 
performed  on  the  unidirectional  specimens  are  planned  for  this  speumen. 


c. 


IlBl  Dependent  Properties  of  Fiber  Composites 
(Dr.  R.  E.  Smith,  Union  Carbide) 


..  e  Orthogonal  ply  graphite -fiber,  epoxy-matrix  composites  are  fre- 

iaVC  Sma,U  CrackS  rUnning  thrOUgh  each  P1/  even  Prior  to  testing.  (2,  34) 
h"  XSe  n  Crt  iS  SeeIi  in  FigUre  41-  -Although  these  cracks  appear  to“  ~ 

have  little  influence  on  the  composite  ultimate  strength  or  Young's  modulus  in  the 

or  strL°nS  °f  \lheT9>  their  effect  on  time  dependent  properties  (fatigue,  creep, 
m on  rt  r«pture),  on  resistance  to  corrosive  environment,  or  on  off-axis  static 
pioperties  has  not  been  studied  in  detail.  Since  plies  with  cracks  are  sometimes 
c  side  red  to  have  already  failed,  the  effective  utilization  of  graphite-fiber  com- 

dep®nf  on  afUth°rOUgh  understanding  of  the  interaction  of  these  cracks 
with  the  stress  state  of  the  composite. 


Figure  41.  Cracks  in  an  "As  Fabricated  '  "Thornel"  50  Plate. 

N-22511 
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Doner  and  Novak  (34)  have  studied  the  crack  network  under  condi¬ 
tions  of  repetitive  loading  of  composites  of  "Thornel"  50  in  various  resin 
matrices.  Their  work  was  limited  to  fewer  than  1.4  x  1  03  cycles  on  mate¬ 
rial  that,  because  of  the  low  bonding  strength  of  resin  to  untreated  fiber,  had 
a  large  number  of  cracks  in  the  as -fabricated  state.  The  goal  of  this  project 
was  to  extend  our  knowledge  of  composite  cracking  by  using  "Thornel"  50S 
yarn  which  had  been  treated  to  provide  strong  fiber-matrix  bonding,  and  by 
extending  the  range  of  test  conditions  to  sustained  loading  of  up  to  1000  hrs 
and  to  fatigue  cycling  for  up  to  lx  1 06  cycles.  It  was  hoped  that  the  better 
bonding  character  of  the  treated  yarn  would  give  an  essentially  crack-free 
composite  with  which  to  study  the  growth  and  multiplication  of  cracks  under 
stress.  It  was  expected  that  degradation  of  the  composite  tensile  strength 
and  initial  modulus  with  repeated  stress  could  be  correlated  with  the  increased 
cracking  in  a  useful  way. 

The  composite  plates  fabricated  for  this  program  were  12"  square 
flat  panels  with  7  plies  of  "Thornel"  50S,  ERLA  4617/MPDA  in  a  C-gO"  orthog¬ 
onal  layup  with  the  0°  direction  corresponding  with  the  fiber  direction  of  four 
of  the  plies.  The  properties  of  the  plates  and  plate  constituents  are  listed  in 
Table  IX.  Two  sample  geometries  were  used,  one  for  sustained  loads  and  one 
for  fatigue  loads.  Samples  for  sustained  loading  were  1  1/2  inches  wide,  necked 
to  1/2  inch,  and  10  1/4  inches  long.  Reinforcing  doubler  plates  cut  from  stock 
sheet  steel  were  bonded  to  the  samples,  and  the  sandwich  was  drilled  for  pin 
loading.  The  samples  for  fatigue  loading  were  9  inches  long  and  1/2  inch  wide 
with  fiber  glass  doublers. 


TABLE  IX 

PHYSICAL  PROPER.  TES  OF  PLATES  AND  CONSTITUENTS 


Property 

Plate 

~r~ 

T~ 

T 

4 

5 

6 

Tensile  Strength  (ksi) 

Yarn 

220 

220 

220 

220 

220 

200 

Resin 

19 

19 

19 

19 

19 

19 

Plate  (0°  direction) 

55 

66 

68 

62 

67 

72 

Initial  Modulus  ( 1 06  psi) 

Yarn 

52 

52 

52 

52 

53 

58 

Resin 

.78 

.  78 

.  78 

.  78 

.  78 

.  78 

Plate  (0°  direction) 

15.  2 

15.  3 

14.  5 

15.  7 

18.8 

17.  3 

Shear  Strength  (ksi) 

Yarn 

8.  2 

8.2 

8.  2 

8.  2 

8.  0 

7.  3 

Resin 

15 

15 

15 

15 

15 

15 

Yarn  Volume  Content  (%) 

52 

56 

52 

56 

60 

57 
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The  sustained  load  tests  were  conducted  on  Satec,  Model  D,  creep 
rupture  testers.  Specimens  had  strain  gages  on  each  side  which  monitored 
e  strain  as  the  test  progressed.  The  specimens  were  loaded  to  different 
stresses  for  1000  hours  with  the  exception  of  three  specimens  from  which 

ou6t  !  z  ^ I"  rCmrd  ^  2;  2°’  and  200  hOUrS‘  FdtiSue  cycl^g  was  carried 
out  on  a  Gilmore  Universal  Testing  machine  which  is  of  a  servo-controlled 

design'  Tests  were  run  in  a  tension-tension  mode  at  6  cycles/ sec 
with  the  minimum  load  approximately  1/10  of  the  maximum  load. 

,  ,  .  ,  For  PurPoses  of  determining  desired  stress  levels  in  sustained 

load  or  in  fatigue  and  as  a  basis  for  comparison  of  the  effects  of  dynam- 
loads,  it  was  necessary  to  know  the  tensile  strength  and  modulus  to  be 
expected  of  each  specimen;  however,  the  limited  number  of  specimens 
obtainable  from  each  plate  precluded  a  statistical  evaluation  of  the  tests  in 
his  program  The  tensile  strength  of  the  plates  used  in  the  sustained  load 
ests  were  estimated  by  stressing  one  specimen  from  each  plate  to  failure 
More  measurements  of  modulus  were  included  to  determine  the  effect  of 
s  ress  on  modulus.  In  the  fatigue  program,  more  specimens  were  available 
and  baseline  strength  and  modulus  were  determined  from  the  average  of  three 
tests.  Optical  counting  of  cracks  was  chosen  to  provide  straightforward  com- 
senteHn  existing  data.  (34)  The  results  of  optical  crack  counting  are  pre¬ 
sented  in  this  report  in  terms  of  the  "crack  density,"  which  is  defined  her^ 

-  H  fnTS?ri°f  cJuackf8  °bserved  in  a  section  of  the  specimen  divided  by  the 
product  of  the  length  of  the  specimen  and  the  number  of  plies  being  examined 
The  crack  density  ivill  be  given  in  units  of  c  racks/ mm/ ply.  So  that  the  maicr 
source  of  error  in  optical  counting  might  be  minimized,  i.  e.  ,  finding  a  rep¬ 
resentative  cross-section,  each  specimen  was  cut  into  several  sections 

TlfsT.Zr  VeaT  end  °f  thC  gage  SCCtion  and  from  th*  center. 

permitted  us  to  de,ermine  the  °f 

.  Experimental  Results.  One  of  the  prime  goals  of  this  program 

°f  Cf.acking  in  an  initially  crack-free  composite  Other 
studies  (34)  had  found  cracks  in  all  composites  of  "Thornel"  50  and  resins 

Calb  76fre  ^  ^  gh  temPeratures-  Ail  of  the  plates  fabricated  by  Union 

Carbide  for  this  program  were  essentially  crack-free,  a  condition  which  is 

f°,r  cross‘Piy  graphite-fiber  composites  and  is  probably  due  both  to 

f  reungth  °f  the  ^,rn  <7  ksi  vs-  4  k®i  ^  Reference  34)  and  to 
,  ®  b.r.lcatlon  technique.  Many  sections  of  the  plates  were  checked  metal- 
gr  phically,  and  some  poros.meiry  work  was  also  carried  out,  but  no  meas- 
“h  f.1®  crackicg  was  found.  The  microstructure  of  the  plates  was  excellent- 

The  dlstrioutl°" ,  wds  uniform,  and  there  were  no  cracks  and  no  voids. 

S*  s trenBthYuTfSirineHrfh harac,te.riz<;d  by  measurements  of  the  ultimate  ten- 
SEMl  nf  fr?  !  S)  nd  the  mldu'us  and  bV  scanning  electron  microscopy 

ablC  IX  Sh°WS  the  UTS  and  E  values  for  the 
six  plates  used.  The  stress-strain  curves  for  all  specimens  were  nonlinear 

n  appeared  to  consist  of  two  linear  regions.  At  low  stresses  (<  10  ksi) 
there  was  a  linear  region  with  modulus  E„  the  initial  modulus;  at  high 
esses  (>20  ksi),  there  was  a  second  linear  region  with  a  modulus  E2>Ej. 
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In  order  to  provide  a  basis  for  comparison  for  long  term  and 
cyclic  loading,  the  crack  density  was  measured  in  several  specimens  after 
they  had  been  taken  to  various  stresses  applied  in  the  0"  direction  and  re¬ 
turned  to  zero  stress.  These  "single  cycle"  crack  densities  are  presented 
in  Table  X  and  Figure  43.  Table  X  shows  that  the  crack  densities  correlate 
well  with  the  single  cycle  "strain"  calculated  from  the  applied  stress  and 
the  average  initial  modulus.  The  "strain"  thus  calculated  is  merely  indic¬ 
ative,  since  it  does  not  take  into  account  the  variability  of  sample  properties 
or  the  nonlinearity  of  the  stress-strain  curve.  There  appears  to  be  a 
threshold  for  crack  creation  in  the  neighborhood  of  a  "strain"  of  .  30  percent. 
Below  this  value,  the  crack  density  is  at  or  below  the  threshold  of  detect¬ 
ability;  above  .  30  percent,  the  crack  density  increases  rapidly  with  increas¬ 
ing  strain.  Table  X  shows  the  cracking  only  in  the  90®  plies;  the  crack 
density  in  the  0“  plies,  where  the  stress  is  not  oriented  so  as  to  open  cracks, 
was  generally  below  the  detectability  threshold  and  is  not  listed.  The  nature 
of  the  cracking  is  shown  in  Figure  42,  which  is  a  cross -section  of  a  sample 
listed  in  Table  X  which  had  seen  a  stress  of  93  percent  UTS  with  a  "strain" 
of  .  37  percent.  The  cracking  induced  by  stress  appears  in  every  way  simi¬ 
lar  to  that  which  occurs  during  fabrication. 


TABLE  X 

CRACKING  DUE  TO  SINGLE  CYCLE  STRESS 


Plate  No. 

Single  Cycle 
Stress 

<%  of  UTS)  (ksi) 

Single  Cycle 
Strain 
(%) 

Crack  Density  in 
90°  plies 
(cracks/mm/  ply) 

4 

58 

36 

.23 

.  02 

93 

58 

.  37 

.  17 

99 

61.  £ 

5  .  39 

.37 

3 

74 

50 

.  35 

.  18 

83 

56 

.  39 

.  30 

5 

80 

54 

.29 

.01 

Figure  43  presents  the  crack  density  results  for  the  sustained 
stress  samples.  The  samples  were  tested  in  two  groups.  The  first  group 
was  started  in  August  1969  and  the  second  group  in  December  1969.  There 
is  a  marked  difference  in  the  crack  density  of  the  two  groups.  The  four  sam¬ 
ples  from  IHe  first  group  show  "normal"  ;rack  densities  compared  with  the 
single  cycle  stress  samples  and  are  presented  in  Figure  43.  The  crack  densi¬ 
ties  are  higher  for  a  given  strain,  which  might  be  interpreted  as  additional 
damage  due  to  the  sustained  loading.  The  trend  of  crack  density  with  strain 
is  the  same.  The  results  from  the  second  group  were  completely  unexpected. 
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Figure  42.  Cracks  in  Plate  No.  4  Due  to  One  Cycle  Stress 

of  93  Percent  UTS. 

tv,  ,  .  .  .  N-22509 

sameCota^  ;nhltler  7”  mUCh  !°°  l0W  in  ComParison  with  samples  from  the 
same  plate  in  the  first  group  and  with  the  single  cycle  stress  samples  and  are 

and  fih°Wn  m  l  lg“rf  4?\  The  strain  during  the  sustained  loading  was  monitored 
verv  lmT  d  fiPld  initiaI _  Permanent  deformation  in  the  first  few  hours  with 

aooLximlid?f0Anriatl0ti!>th/  reafter’  The  Permanent  deformation  was  small, 
approximately  60p  inches/ inch  at  stress  levels  of  50  ksi  (75  percent  UTS),  and 

n:‘u^;rrr0n  amORg  grouPs  of  ^P1®8-  The  Change  in  modulus 
due  to  sustained  load  was  measured  and  found  to  be  insignificant. 

.  .  The  results  of  fatigue  cycling  for  1  x  105  cycles  at  various  stresses 

are  also  presented  in  Figure  43  and  the  results  of  cycling  for  a  range  of  number 

rncreCasee8dastres!tant-StreSS  7*  preusented  in  Table  XI-  There  is  a  fendency  for 
dens i  °r  increased  number  of  cycles  to  produce  increased  crack 

density,  but  there  are  two  anomalous  features  which  limit  the  validity  of  anv 
nclusions.  First,  some  samples  show  unusually  low  crack  counts  for  tbe^ 
thereS;oanf  c/cl«  condltions »  Just  as  in  the  sustained  load  samples.  Second, 

sheeted r®quentiy  a  lar8e  discrepancy  between  the  crack  densities  of  samples 
subjected  to  essentially  identical  stress -cycle  conditions. 
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Crack  Density  (cracks/ mm/ ply) 


Figure  43.  Crack  Density  Produced  by  Various  Strain-Time  Conditions. 

G-700137 
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TABLE  XI 

CRACKING  DUE  TO  FATIGUE  CYCLING  AT  CONSTANT  STRESS 


Maximum  Stress  Maximum  Strain  Number 
Specimen  (ksi)  (%)  of  Cycles 


5-13 

54 

5-19 

54 

5-16 

54 

5-8 

54 

5-7 

54 

5-11 

54 

6-2 

57.  5 

6-14 

57.5 

6-1 

57.  5 

6-7 

57.  5 

6-15 

57.  5 

6-9 

57.  5 

6-10 

57.  5 

6-  12 

57,  5 

6-4 

57.  5 

6-16 

57.  5 

.29 

1 

.  29 

10 

.  29 

1  x  103 

.  29 

1  x  104 

.  29 

1  x  10 

.  29 

i  x  106 

.  33 

10 

.  33 

10 

.  33 

1  x  103 

.  33 

1  x  10J 

.  33 

1  x  104 

.  33 

I  x  104 

.  33 

1  x  10s 

.  33 

i  x  io3 

.  33 

lx  10* 

.  33 

lx  10‘ 

Crack  Density 
(cracks/ mm/ ply) 

.  01 
.  11 
0.  0 
.  04 
.51 
.67 
.03 
.03 
.  12 
.  39 
.25 
.20 
.  55 
.83 
1.65 
.69 


,The  tensile  strength  and  modulus  are  not  much  affected  by  fatieue 
cycling,  as  shown  in  Table  XII.  There  may  be  a  modest  decrease  in  UTS  for 
the  high  strain  cycling  and  an  increase  for  the  lower  strain  conditions.  The 
small  number  of  samples  tested  makes  the  apparent  increase  in  UTS  for 
.?**  5  aftfr  fLati6ue  statistically  inconclusive,  but  the  result  is  consistent 
with  that  of  other  investigators.  (35)  In  any  case,  fatigue  cycling  does  not 
adversely  affect  the  tensile  strength  of  the  composite 


TABLE  XII 

MODULUS  AND  STRENGTH  OF  FATIGUED  SPECIMENS 


Plate  No. 

Number  of 
Cycles 

Maximum 

Stress 

(ksi) 

Maximum  Tensile 
Strain  Strength 
(%)  (ksi) 

Initial 
Modulus 
(10*  psi) 

3 

baseline 

68 

14  5 

fatigue 

C 

1  x  10s 

56 

.  39 

66 

i5;  3 

D 

baseline 

...  - 

••  as  _ 

65, 

67, 

69 

18.  8 

fatigue 

1  x  104 

58 

.31 

69, 

71, 

74 

19.  0 

1  x  105 

58 

.  31 

66, 

68, 

71 

18.  2 

1  x  106 

58 

.  31 

72, 

75, 

76 

19.  3 
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an  attempt  to  evaluate  the  crack  density  of  samples  which  had 
tailed,  we  found  that  fractured  specimens  had  a  very  nonuniform  crack  count. 

I  he  crack  density  was  very  high  near  the  point  of  fracture  and  fell  off  steadily 
as  a  function  of  distance  from  the  fracture.  Intact  specimens  from  the  sus¬ 
tained  loading  and  fatigue  cycling  occasionally  also  showed  striking  honuni- 
formities  of  crack  density,  varying  by  a  factor  of  2,  but  there  was  no  consis¬ 
tent  trend  for  cracking  to  be  high  near  the  grips,  for  instance. 

rhe  problems  of  abnormally  low  crack  densities,  lack  of  repro¬ 
ducibility  of  crack  densities,  and  nonuniformity  of  crack  densities  make  dif¬ 
ficult  the  meaningful  interpretation  of  data  generated  during  this  program. 

These  three  problems  are  very  likely  due  to  one  or  more  of  the  following' 
causes:  intiinsic  sample  variability,  sample  preparation,  variability  in  test 
conditions,  or  variability  in  post-test  examinations.  No  single  explanation 
seems  entirely  sufficient.  Testing  conditions  were  monitored  as  carefully 
as  possible.  Sample  preparation  and  examination,  which  might  create  cracks, 
would  hardly  generate  crack  densities  which  vary  smoothly  with  test  strain. 
Sample  variability  is  inadequate  to  explain  the  low  cracking  in  all  of  the  second 
group  of  sustained  load  specimens. 

The  lack  of  reproducibility  of  crack  densities  in  samples  fatigued 
under  the  same  nominal  conditions  may  be  easier  to  explain.  Although  the 
samples  were  loaded  to  a  constant  maximum  stress,  Table  XII  indicates  that 
a  constant  maximum  strain  is  more  relevant  to  cracking.  The  rapid  variation 
of  crack  density  with  strain  coupled  with  the  normal  sample-to-sample  varia¬ 
tion  in  modulus  (and,  hence,  strain  for  a  given  stress)  could  well  account  for 
the  lack  of  reproducibility  of  crack  density.  A  crude  estimate  from  Figure  43 
indicates  that  a  variation  in  strain  of  10  percent  is  capable  of  producing  a  varia¬ 
tion  of  a  factor  of  two  in  the  crack  density  for  single  cycle  stress  and,  pre¬ 
sumably,  also  in  the  fatigue -induced  crack  density. 

^  ‘e  n°nuniformity  of  crack  densities  in  fractured  samples  can  be 
explained  in  two  ways.  In  fracture,  the  elastic  strain  energy  is  released 
suddenly  as  a  shock  wave  which  propagates  down  the  sample,  gradually  dim¬ 
inishing  in  intensity.  The  cracking  would  be  indicative  of  the  shock  amplitude. 
The  principal  problem  with  this  interpretation  is  that  the  crack  densities  far 
from  the  fracture  were  observed  in  most  cases  to  be  lower  than  those  expected 
from.  Table  X.  This  conditiomwould  seem  to  indicate  that  the  crack  density 
can  build  up  locally,  around  a  stress  concentrating  defect,  for  instance,  and 
relieve  the  strain  in  the  90°  layers  elsewhere  in  the  sample.  The  nonuniform¬ 
ities  of  the  crack  density  of  intact  fatigue  specimens  may  be  related  to  this 
effect  and  may  indicate  the  probable  region  of  ultimate  failure  of  the  specimen. 

In  view  o.f  the  crack  anomalies,  there  is  little  in  the  way  of  analysis 
that  can  be  done  on  the  cracking  due  to  sustained  loads.  Samples  may  or  may 
not  have  cracking  for  reasons  which  may  only  be  guessed..  .  The  primary 
effect  of  the  sustained  stress  seems  to  be  to  lower  the  apparent  threshold  strain 
for  crack  creation  from  .  3  percent  for  single  cycle  stress  to  less  than  .  2  percent 
for  sustained  stress.  Since  cracking  varies  with  strain  rather  than  stress,  it 
is  not  possible  to  specify  the  crack'density  for  a  plate  at  a  given  stress,  such 
as  the  design  limit  of  2/3  UTS;  however,  for  a  plate  typical  of  this  program 
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(6b  ksi  T  rS  and  18  x  106  psi  modulus),  the  crack  density  after  1000  hours  at 
2/3,UTS  would  be  approximately  .  30  cracks/ mm/ ply  (or  fewer,  if  the  unknown 
conditions  causing  low  crack  densities  are  met).  The  permanent  deformation 
in  these  composite  specimens  is  small  and  correlates  reasonably  well  with 
either  strain  or  stress.  In  particular,  there  are  no  obvious  anomalies  in  the 
second  group.  This  fact  would  appear  to  indicate  that  the  permanent  deform¬ 
ation  is  not  related  to  the  creation  of  cracks. 

Fatigue  cycling  does  not  appear  to  have  a  strong  effect  on  the 
threshold  strain  for  crack  creation  but  does  strongly  affect  crack  density 
above  the  threshold.  Thus,  for  the  "typical"  specimen  discussed  above, 
fatigue  for  10  cycles  at  2/3  UTS  would  create  few  cracks.  This  result  was 
inadvertently  verified  early  in  the  program  when  the  fatigue  machine  failed 
to  follow  the  set  load  and  a  run  was  made  at  a  stress  <2/3  UTS.  No  cracks 
were  found  in  this  specimen. 

The  trend  in  crack  density  in  this  program  may  be  discussed  in 
comparison  with  the  results  of  Doner  and  Novak.  (34)  They  find  that  crack 
density  increases  with  stress  but  that  the  density  is  not  the  same  for  different 
plates  for  the  same  stress  as  a  percentage  of  UTS.  This  result  is  consistent 
with  the  finding  of  this  program  that  cracking  is  related  to  strain.  Doner  and 
Novak  also  find  that  cracking  occurs  mainly  during  the  first  cycle,  whereas 
this  study  shows  that  cracking  increases  steadily  with  number  of  cycles.  This 
disagreement  is  probably  due  largely  to  the  different  materials  involved.  In 
this  study,  the  crack  density  at  106  cycles  and  a  stress  of  80  percent  UTS  was 
less  than  that  which  Doner  and  Novak  found  after  one  cycle  at  40  percent  UTS 
in  their  lower  shear  strength  material.  In  any  case,  cyclic  loads  at  stresses 
up  to  85  percent  UTS  and  for  up  to  106  cycles  appear  to  produce  no  degradation 
of  tensile  properties  (Table  XII  and  references  34  and  35).  In  fact,  there  may 
even  be  a  slight  strengthening  effect  of  fatigue,  perhaps  due  to  the  relief  of 
thermal  stresses. 

The  crack  inducing  character  of  fatigue  differs  from  that  of  sus¬ 
tained  stress  (Figure  43).  Within  the  limits  already  discussed,  the  sustained 
stress  appears  to  lower  the  threshold  strain  for  crack  creation  but  does  not 
affect  the  slope  of  the  creek  density-v.s-test  strain  curve  much.  Fatigue  for 
10  cycles,  however,  does  not  appear  to  affect  the  crack  threshold,  but  it 
does  appear  to  increase  the  slope  of  the  cracking-strain  curve  dramatically. 

D.  Notch  Sensitivity  of  Fiber  Composites 

(Dr.  IT!  Et  Smith,  Union  Carbide) 

The  objective  of  this  work  was  a  comparative  evaluation  of  the 
nolch  sensitivity  of  "Thornel"  50S,  ERLA  4617  epoxy  and  P-1700  polysul'one 
composites.  Since  polysulfone  resins  are  tougher  and  have  much  higher 
elongati-ow  to  failure  than  the  relatively  brittle  epoxies,  it  was  expected  that 
the  polysulfone  composites  would  be  less  notch  sensitive.  Polysulfone  compos¬ 
ite  plates  became  available  only  at  the  very  end  of  this  program;  thus,  only 
a  small  number  of  experiments  were  performed. 
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A  12-inch  square  7-ply  (0°,  90°)  plate  of  "Thornel"  5  OS/ ERL  A 
4&17/MPDA  composite  was  cut  into  seven  specimens  according  to  standard 
specifications  (MIL -Handbook-  1 7,  pg.  2-29).  T’.e  specimens  were  1  1/2  inch 
wide,  necked  to  1/2  inch  by  arcs  of  4  inch  radium.  A  1/8  inch  hole  was  drilled 
into  the  center  of  the  necked  section  in  four  of  the  specimens.  The  failure 
stress  was  calculated  from  the  failure  load,  and  the  remaining  cross-sectional 
area  was  determined  from  a  micrometer  measurement  of  neck  width,  thickness 
and  hole  diameter.  Reinforcing  doublers  similar  to  those  used  in  the  fatigue- 
cracking  program  were  employed.  In  addition,  a  12-inch  square  7-ply  (±45°) 
plate  was  prepared  from  the  same  prepreg  and  cut  similarly  to  test  the  varia¬ 
tion  of  notch  sensitivity  with  ply  orientation.  The  fabrication  of  polysulfone 
plates  is  reported  in  Volume  I  (Materials  Research)  of  this  Report.  Since  it 
was  difficult  to  prepare  large  polysulfone  plates,  it  was  not  possible  to  cut 
standard  specimens.  We  decided  to  try  samples  8  inches  long  by  1/2  inch  wide, 
®  the  fatigue  samples,  with  the  same  doublers  as  in  the  fatigue  pro¬ 

gram,  One  plate  (PS-08270)  of  "Thornel"  50S/ polysulfone  was  chosen  for 
this  work,  based  on  this  microstructural  examination  (Figure  44).  Six  speci¬ 
mens  were  cut  from  this  12-inches  by  4-inches  plate.  Three  of  these  samples 
had  holes  drilled  into  the  center  of  the  gage  section.  So  that  this  change  in 
specimen  geometry  cv  uld  be  evaluated,  three  "Thornel"  50S,  ERLA4617  speci¬ 
mens  of  this  straight-sided  geometry  were  also  tested  for  notch  sensitivity. 
These  samples  were  cut  from  fatigue  Plate  No.  6.  (Table  IX,  Section  VI-C  of 
this  report). 


Figure  44.  Microstructure  of  Plate  PS-08270. 
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presents  the  results  of  the  notch  sensitivity  testing. 
The  necked  epoxy  composite  specimens  lost  30  percent  of  their  strength 
when  tested  along  a  ply  direction  and  less  than  2  percent  of  their  strength 
when  tested  at  45c  to  the  ply  directions.  The  straight-sided  geometry  epoxy 
composite  specimens  tested  along  a  ply  direction  lost  18  percent  of  their 
strength.  All  of  the  epoxy  specimens  failed  in  a  "normal"  manner:  the  0° 
failures  were  brittle  with  irregular  failure  surface  and  the  45*  failures 
failed  along  the  fibers  leaving  a  sawtooth  surface.  The  polysulfone  com¬ 
posites  appear  to  have  no  notch  sensitivity  when  tested  along  a  ply  direc¬ 
tion,  a  conclusion  which  is  somewhat  tentative  because  the  unnotched  (base¬ 
line)  specimens  all  failed  in  the  grips;  thus,  the  reported  baseline  failure 
stress  is  only  a  lower  limit.  However,  the  failure  stress  of  64  ksi  of  the 
notched  specimens  is  close  to  that  expected  on  a  rule-of-mixture  basis. 


TABLE  XIII 


COMPOSITE  NOTCH  SENSITIVITY 


Control 

Notched 

Specimen 

Strength 

(ksi)  No.  of  Tests 

Strength  No.  of 
(kai)  Tests 

Strength 
Loss  (%) 

Epoxy 


Necked  (0°,  90°) 

51 

3 

35.2 

4 

31 

Necked  (±45°) 

14.6 

3 

14.  4 

4 

2 

Straight  (0° ,  90°) 
Polysulfone 

71.9 

3 

59.  1 

3 

18 

Straight  (0° ,  90°) 

61.2 

3 

64.  3 

3 

0 

E. 


Influence  of  Defects  on  Composite  Properties 
(Dr.  T~.  Weng,  Union  Carbide)  C 


Earlier  in  this  program  (3),  four  1 2 -inch  square  panels  which 
contained  a  variety  of  known  intentional  defects  were  fabricated  for  use  in 
evaluating  nondestructive  test  methods.  The  panels  were  made  from 
Thornel"  40/ERL  2256/MPDA  and  had  a  balanced  orthogonal  lay-up. 

Two  panels  were  made  with  four  plies  and  two  panels  were  made  with*  12 
plies.  Defects  consisted  of  Mylar  and  off-axis  yarn  inclusions  between 
plies,  butt  joints,  lap  joints,  and  resin-rich  and  resin-starved  areas. 

These  defects  were  thought  to  be  representative  of  those  which  might  occur 
in  composite;  fabrication. 


After  nondestructive  testing  was  carried  out  at  AVCO  Corporation  (36) 
test  specimens  were  cut  from  these  plates  to  evaluate  the  effects  of  these  detects"  * 
on  the  tensile,  flexural,  and  compression  strength  of  the  composite.  The  sonic 
modulus  of  the  specimens  was  also  measured  prior  to  the  static  destruction  tests. 
The  test  results  are  listed  in  Table  X7.V. 

Butt  joints  and  lap  joints  reduced  the  tensile  strength  by  approximately 
25  percent.  Double  yarn  inclusion  had  no  apparent  effect  on  tensile  strength. 

Butt  joints  and  lap  joints  appear  to  have  only  a  minor  degrading  effect  on  flexural 
strength,  but  caution  may  still  be  required  because  the  effects  on  flexural  strength 
depend  not  only  on  the  type  of  defect  but  aloo  on  its  location.  The  defects  in  the 
flexural  specimens  were  located  at  approximately  one-fourth  of  the  thickness  from 
a  surface,  but  it  was  not  known  whether  they  were  located  in  the  tencile  or  com¬ 
pressive  region  of  the  sp^eimens.  The  effect  of  defects  on  compressive  strength 
was  smaller  than  that  on  nsile  strength.  It  is  probable  that  stress  concentra¬ 
tion  from  the  defects  are  smaller  under  compressive  lc'ads  th&u  under  tensile 
loads. 


TABLE  XIV 

EFFECT  OF  KNOWN  DEFECT  ON  COMPOSITE  STRENGTH 


Plate 

Sonic 

Strength 

Modulus 

Strength 

Retention 

Numb 

Test 

Defect 

(10‘  psi) 

(psi) 

1%) 

Control  (2)  {a) 

11.6 

64,200 

100.0 

3 

Tension 

Butt  Joint  (2) 

11.6 

46,700 

72.  7 

4  ply  (0° , 90°) 

Double  Inclusion  Yarn(l)  11.8 

61,300 

95.4 

2 

Control  ( 1 ) 

11.4 

59,000 

100.  0 

12  ply  (0°  ,  90°) 

Tension 

Lap  Joint  (2) 

11.0 

46,200 

78.  3 

1 

Flexural ^ 

Control  ( 1 ) 

9.  5 

67,500 

100.  0 

12  ply  (0° , 90°) 

Lap  Joint  (2) 

9.  2 

65,200 

96.5 

Flexural 

Control  (1) 

11.6 

78,400 

100.0 

2 

Resin  Starved  (1) 

12.  1 

83, 100 

105.9 

12  ply  (0°  ,  90° ) 

Butt  Joint  ( 1 ) 

10.  7 

72,400 

92.  3 

Control  ( 1 ) 

10.  7 

42,400 

100.0 

4 

4  ply  (0°  ,  90°) 

Compression 

(c)  Butt  Joint  ',) 
Lap  Joint  (2) 

10.  3 

12.  3 

39,200 

37,400 

36,300 

92.4 
86.  2 

Mylar  Inclusion  (1) 

9,6 

85.6 

jfjNumber  in  parentheses  is  number  of  test  sample. 
.  (Single-Point  loading. 

'C  Fed.  Test  Method  Std.  No.  406. 
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~'raCt^f.e  Stren?th  °f  Composite  Cylinders  Under  Combined  Stresses 
(Dr.  X.  Weng,  Union  C arbide )  “  “ - 

,  Experimental  data  on  the  stress -strain  properties  and  fracture 
of  composites  under  combined  stress  loading  are  needed  for 
evaluating  the  accuracy  of  various  theories  of  fracture  for  fiber -reinforced 

composites  *  ^  ™*6ily  avaiIable  to  date  on  graphite  fiber 

f  Therefore,  an  experimental  investigation  was  conducted  to 
generate  fracture  data  under  combined  stresses  with  which  the  various 
theories  of  fracture  may  be  compared. 

stresses  an S,tTe^th  under  biaxial,  torsional,  and  combined 
stresses  and  the  stress -strain  properties  of  "Thornel"  50S/EREB  4617/MPDA 

no  hre  ^een  detelrmined'  The  composites  were  of  four-ply  (90%  0° 

0  ,  90  )  construction,  with  a  fiber  content  of  55  to  58  volume  percent. 

.  Torsional  and  combined  stress  tests  were  conducted  on  cylindrical 

renorwS?  TV,  i  C°m°ined  stress  test  apparatus  described  in  a  previous 
90°  “f  lay;uP  Pattern  of  the  cylindrical  specimens  was  such  that  the 
LJ  f\ltl0n  Coincided  wlth  the  tangential  direction  of  the  cylinder.  The 
ends  of  the  specimens  were  reinforced  with  liber  glass  prepreg  cloth  and 
were  attached  to  the  steel  fixtures  with  epoxy  adhesive  a*  shown  in  Figure  45 
The  dimensions  of  the  specimens  were:  inside  diameter,  2.  01  in.  ;  outside 
lameter,  2.  08  in.  ;  gauge  length,  2.  50  in.  ;  and  total  length,  7.  0  in. 

fL  c_  Ten  cylindrical  specimens  were  fabricated  with  material  from 

the  same  prepreg  run.  Two  specimens  were  tested  in  pure  torsion  and 

selected"^  teStedvUndHr  COmbined  ladings.  The  load  conditions  were 

*  Permit  a  determination  of  the  approximate  shape  of  the  fracture 

rlsuHs  of  th  encomP.assing  the  shear  -tensile-tensile  stress  states.  The 

Ire  g^en  in  Tlb^Xv"  "  ^  StreSS  t68t8  °"  *e  ^lindrical  ^ecimens 

„ _ _  .  In  Edition  to  the  combined  stress  tests,  the  uniaxial  tensile  and 

0°  °°  direct.ion  were  measured  with  (90°,  0°, 

,  vu  )  and  (0  ,  90  ,  90  ,  0  )  composite  plates,  and  the  short  beam  shear 

The  r^dtro?SVJirtd  Tith  a  f°Urteen.-P1y  unidirectional  composite  plate, 
l  e  results  of  these  tests  are  summarized  in  Table  XVI.  Strength  data 

E£J».  Ta“es  XV  except  the  short  beam  shear  Strang*,  are 

plotted  in  a  three-dimensional  graph  with  coordinates  of  shear,  tangential 

denoted8wUhthl  tha^T  7"  in/igur*  46'  The  uniaxial  strength  data  kre 
denoted  with  the  shaded  circles  and  combined  strength  data  are  denoted  with 

t  KPuin  ':ircles'  Curves  shown  in  Figure  46  were  drawn  to  indicate  the 
probable  fracture  strength  surface  of  the  composite. 

?ne  C7ndrical  sPecimen  (09140-1 )  was  tested  under  internal 
pressure  and  axial  compressive  load.  The  objectives  of  this  test  were  to 
simulate  the  condition  of  a  hoop  tensile  test,  where  the  magnitudes  of  the 

lenVil'}' .T  “d  r/dial*«esses  very  small  in  compariCm  the  ta„! 
gential  stress  and  to  examine  the  effects  of  end  constraints  on  the  specimen. 


-  107  - 


External  End  Holder 
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Figure  45.  Combined  Stress  Test  Specimen  Assembly. 


r,  Shear  Stress 


o"2»  Longitudinal 
Stress 


Stress  Unit:  103  psi 


Figure  46.  Fracture  Strength  Surface  for  Four-Ply  (90°  06  0°  90°) 

"Thornel"  50S/ERLB  4617  Composite.  ’ 

G-710026 
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^Tested  in  the  0*  direction 


The  specimen  was  first  subjected  to  an  internal  pressure  of  480  psi  prior  to 
the  simulation  of  the  hoop  tensile  condition.  Tangential  and  longitudinal 
strains  at  the  center  of  the  gauge  section  and  at  the  edge  of  the  fiber  glass 
reinforcement  were  measured.  The  magnitude  of  the  tangential  strain  at 
the  end  was  about  95  percent  of  that  at  the  center.  At  fracture  under  internal 
pressure  and  axial  compressive  lr,ad.  the  tangential  stress  was  47,  700  psi 
and  the  longitudinal  stress  was  1100  p3i;  therefore,  the  condition  of  the  hoop 
tension  was  approximately  obtained.  Tangential  strains  at  the  center  and 
end  of  the  gauge  section  were  found  to  be  approximately  the  same  at  fracture. 
Therefore,  the  existing  end  constraints  appeared  to  have  no  significant  effect 
on  the  fracture  strength  of  the  composite  since  the  tangential  stress  of 
47,  700  psi  wa  '  about  15  percent  greater  than  the  uniaxial  tensile  strength 
determined  in  the  0°  direction  with  the  (0°,  90°,  90°,  0°)  plate. 

The  strengths  of  the  pure  torsion  specimens  were  low  in  compar¬ 
ison  with  the  remaining  combined  stress  test  data.  These  specimens  had 
some  wrinkles  on  the  outside  surface  in  the  gauge  section  and  these  defects 
could  have  affected  their  strength.  A  more  realistic  value  of  the  shear 
strength  has  been  estimated  for  the  cylindrical  specimens  by  extrapolation  of 
the  fracture  strength  of  the  cylindrical  specimens  tested  under  combined 
stress  loadings  to  the  condition  of  pure  torsion  (to  the  limit  of  zero  axial 
load  and  zero  internal  pressure).  The  extrapolated  fracture  surface  is  shown 
m  Figure  46,  and  the  estimated  shear  strength  was  7600  psi.  However,  this 
ef .valu.e  of  shear  Strength  is  still  lower  than  the  shear  strength  value 
of  9200  psi  which  was  measured  by  the  short  beam  shear  test  on  the  fourteen- 
p  y  unidirectional  plate.  It  is  not  known  whether  these  two  shear  strength 
values  should  be  identical,  but  this  discrepancy  could  be  due  partially  to  the 
difference  in  the  qualities  of  the  cylindrical  specimens  and  the  plate.  Micro¬ 
photographs  of  the  pure  torsion  cylindrical  specimens,  a  typical  cylindrical 
specimen  subjected  to  combined  stresses,  and  the  fourteen-ply  unidirectional 
plate  are  shown  in  Figures  47  to  50.  Figure  47  is  a  microphotograph  of 
specimen  08100-1  taken  along  the  tangential  direction  of  the  cylinder. 

Figure  48  is  a  mjcrophotograph  of  specimen  08120-1  taken  along  the  longitu- 
mal  direction  of  the  cylinder.  Figure  49  shows  a  microphoto  graph  of 
specimen  08260-2  taken  along  the  tangential  axis  of  the  cylinder;  this  micro - 
photograph  illustrates  the  quality  of  the  typical  combined  stress  test  cylinders. 
Figure  50  is  a  microphotograph  of  the  fourteen-ply  unidirectional  plate.  These 
microphotographs  indicate  that  the  quality  of  the  cylindrical  specimens  has 
been  improved  markedly  from  that  described  previously.  (5)  In  the  pure 
torsion  cylinders,  there  were  more  resin  rich  areas,  voids  between  plies, 
and  cracks  within  plies  than  desired.  The  typical  combined  stress  test 
cylinders  contained  only  a  small  amount  of  these  defects.  The  quality  of  the 
fourteen-ply  unidirectional  plate  was  excellent.  It  is  known  that  there  are 
inconsistencies  among  the  shear  strength  values  measured  by  short  beam 
shear  test  and  pure  torsion  tests  on  tubular  specimens  and  small  solid  rods, 
and  this  problem  should  be  investigated  further. 
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Figure  47.  Microphotograph  of  Pure  Torsion  Test  Cylinder 
(08100-1),  Tangential  View.  (5 OX  Magnification) 


Figure  48.  Microphotograph  of  Pure  Torsion  Test  Cylinder 

(08120-1),  Longitudinal  View.  (5 OX  Magnification) 


Figure  49.  Microphotograph  of  Combined  Stress  Test  Cylinder 
(08260-2),  Tangential  View.  (50X  Magnification) 


Figure  50.  Microphotograph  of  14(0°)  Unidirectional 
Plate  (T50S-091 00-1).  (30X  Magnification) 
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The  preliminary  fracture  surface  shown  in  Figure  46  in  -tee 
that  the  ability  of  the  composite  cylinders  to  withstand  torsional  loads 
simultaneously  with  internal  pressure  and  axial  tension  did  not  decrease 
appreciably  until  the  magnitudes  of  the  internal  pressure  and  axial  tension 
became  relatively  large.  This  characteristic  of  the  composite  is  believed 
to  be  an  inherent  characteristic  of  the  lay-up  pattern. 

Photographs  of  typical  fractured  specimens  are  shown  in  Figure 
51.  Shear  failure  in  and  between  plies  often  occurred  under  pure  torsional 
load.  A  longitudinal  crack  occurred,  as  expected,  when  a  specimen  was 
subjected  to  internal  pressure.  A  fracture  mode  similar  to  that  under  pure 
torsional  load  was  obtained  for  a  specimen  tested  under  torque  and  internal 
pressure,  because  the  specimen  was  tested  under  a  condition  in  which  torque 
was  dominant.  The  fracture  mode  for  specimens  tested  under  combinations 
of  internal  pressure,  torque,  and  axial  load  differed  and  depended  upon  which 
loading  was  dominant  in  the  test. 

In  future  experiments,  the  fracture  strength  in  other  stress 
states,  the  effects  of  loading  paths  on  fracture  strength,  and  other  lay-up 
patterns  should  be  investigated.  It  was  not  possible  within  the  remaining 
time  of  the  contract  to  compare  the  fracture  surface  determined  in  the 
present  experiments  with  fracture  surfaces  predicted  analytically  by  current 
theories  of  failure,  but  the  data  should  be  of  interest  for  such  comparisons. 
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SECTION  VII 


REINFORCED  CUTOUTS  IN  COMPOSITE  STRUCTURES 
L.H.  Kocher,  Bell  Aerospace 


virtuanyUtall°fairf™*S  SiZSS  and  shapes  are  required  in 

lf”ssrr  r  'rv1™11' ,s°s2s*i- 

desianino  *  equipment  and  mechanisms.  Analysis  methods  for 

illicit  a «el?fS7wSu  establishedt0UpS  ln(”4allfc  i-M 

structures,  howeve^  dlsi^lSdeXnes 
Af  Objective  and  Approach 

correlate0 tech StUfY  uere  to  develop  and  experimentally 
around  un£2info?ced  L StfeSS  and  strain  distributions 
The  prinarrana?ytioal  tooi  :af  Ben's  d^retelllmonr  Sh?PeS‘ 

cutout  size*?  anri  thsnne  9  mu  n  generalized  requirements  for 
of  test  specimens:  P  *  Th®  investl9ati°n  considered  two  types 


(1) 


(2) 


pla^es  Vith  circular  cutouts  suojected 
to  uniaxial  tensile  loading,  and 

A  hat-section  stringer-stiffened  panel  with  a 
loading. CUt°Ut  subjected  to  uniaxial  compressive 


B. 


Fabrication 

Initial  Studies 


specimens  were  fabricated  ^with  r6~?ly  unidirectional  plate 
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Corpora,' ion  on  the  prepreg 


characteristics : 

Fiber  56.8 

Resin  41.3 

Finish  0.7 

Volatiles  1.2 


100.0 


indicated  the  following  average 

percent  by  weight 

percent  by  weight 
with  no  volatiles 

percent  by  weight 

percent  by  weight 


Yarn  strand  strength 
Yarn  strand  modulus 
Yarn  density 
Yarn  area/end 
Resin  density 
No.  of  ends/inch 


247  x  10  psi 
6 

51.6  x  10  psi 

1.724  g/cc 

•  "  -5  2 

6.24  x  10  in. 

1.27  g/cc  =  0.0458  lb/ in, 

72 


Torsion  rod  shear 
strength 


6.3  x  10  psi 


Measurements  at  Bell  on  the  prepreg  indicated  43.5  percent  resin 
by  weight  and  volatiles  averaged  3.3  percent.  The  difference 
between  UCC  and  Bell  volatile  contents  was  subsequently  traced  to 

tnno  f^eren?e  in.test  techniques.  (Union  Carbide:  15  minutes  at 
'00oc  in  still  air  oven;  Bell:  15  minutes  in  still  air  oven  at 
160  C.  )  The  4  x  6  in.  panels  were  subsequently  tested  for  in¬ 
plane  notch  sh«iar  and  transverse  tensile  strengths  (37).  Table  XVII 
summarizes  the  test  data,  Table  XVIII  defines  the  bleeder  systems 
investigated.  In  all  cases  Union  Carbide  Corporation's 
recommended  cure  schedule  was  used:  * 


Beat  up  to  20C°F  in  not  less  than  15  minutes. 

Fold  200-210  F  for  60  minutes  under  pressure. 

Increase  temperature  to  240° F  in  not  less  than  15 
minutes,  hold  240-250  F  for  90  minutes. 
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TEST  RESULTS  OF  * ‘THORNEL*  *  50S  UNIDIRECTIONAL  LAMINATES 
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Increase  temperature  to  330^ F  in  not  less  than  30 

minutes,  hold  330-340° F  for  90  minutes. 

Cool  to  ambient  under  pressure  and  vacuum. 

The  cure  pressures  were  nominally  400  psi  for  Panels  1  and  2  and 
100  psi  for  Panels  3-7. 

Although  there  are  not  sufficient  data  on  which  to  base 
firm  comparisons,  it  appeared  in  Table  XVII  that  the  highest 
average  shear  strength  (5.22  ksi)  was  obtained  with  the  use  of 
the  No.  3  bleeder  system,  i.e.(,  with  the  rubber  dam/nylon  bleeder 
materials.  Additionally,  the  No.  3  system  resulted  in  the  best 
laminate  surface  finish.  Higher  transverse  tensile  strengths 
(3,95  and  4.25  ksi)  were  obtained  with  the  No.  2  and  7  systems. 
The  latter  transverse  tensile  strength  increases  (compared  with 
the  3. 1-3.4  ksi  levels  with  the  other  systems)  may  cr  may  not  be 
as  significant  as  the  shear  strength  differences.  In  any  event, 
since  composite  shear  strength  and  surface  finish  are  important 
in  cutout  studies,  the  No.  3  system  was  selected  for  test 
specimen  fabrication. 

2 .  Test  Specimens 

All  test  specimen  plates  were  of  a  fivaiply  (±15,  90, 

T15°)  construction,  selected  for  reasons  of  balance  (cured  flat 
laminates)  and  simulation  of  the  fuselage  component  skin. 
Stringers  for  the  stiffened  panels  employed  a  four-ply  (+10, 
?1(?)  construction,  with  the  same  cross-section  as  that  used  for 
the  fuselage  component.  "Thornel"  50S/4617  prepreg  was  used  for 
all  specimens.  The  unstiffened  tensile  panels  were  6  x  12  in. 
plates  with  measured  VFs  of  55  to  57)5.  The  compression  pa.ial  was 
14.6  in.  wide  x  24  in.  long  with  a  skin  VF  of  55%  and  average 
stringer  VF  was  60)5.  Specimens  were  fabricated  as  follows. 

Unstiffened  Tension  Panels  with  Circular  Cutouts 

(a)  Control  -  Bonded  to  loading  plates  and  tested  as 
shown  in  Figure  52(a). 

(b)  Unreinforced  Cutout  -  A  one-inch  diameter  circular 
cutout  was  machined  in  center  of  panel  (See  Figure 
52(b)). 

(c)  Externally  Bonded  Steel  Reinforcement  -  A  tapered 
4130  steel  disk  having  a  3  in.  O.D.  and  1  in.  I.D. 
was  machined  and  bonded  to  one  side  of  the  test 
panel  by  using  Epon  921  adhesive  and  a  room 
temperature  cure.  Bond  line  thickness  was  10 
mils.  (See  Figure  53.  ) 
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b.  Unreinforced  Cutout 

Figure  52.  Control  and  Unreinforced  Cutout  Test  Panels. 


Figure  53.  Externally  Bonded  Steel 
Reinforced  Cutout  Panel. 
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(<i)  Integral  Composite,  Specimen  (1)  -  Ten  additional 
p  ies  of  qraphite  prepreq  were  interspersed 
between  the  5  plies  of  the  basic  panel  in  the 
region  of  the  cutout.  Each  plv  of  the 

reinforcement  was  trimmed  to  a  different  size  to 
achieve  a  tapering  effect  with  the  thickest  point 
nearest  the  edge  of  the  cutout.  The  resulting  15- 
ply  orientation  at  the  cutout  was  (±15.  ±30  0 

2?/  »  ^ #  '  0  •  90 ,  0,  f  30,3-1 5° ) .  See  Figure 

54(a). 

(e)  Integral  Composite  Specimen  (2)  -  The  same 

fabrication  procedure  and  reinforcement  plv 
orientations  ware  used  as  in  the  integral 

Composite  Specimen  (1)  design,  except  that  the 
configuration  was  changed  from  a  souare  patch  to  a 

butterfly"  type.  See  Figure  54(b). 

(f)  Integral  Titanium  and  Stainless  Steel  Shims  -  Four 

shims  of  titanium  and  stainless  steel  were 

inserted  between  the  basic  5-plv  layup  in  the 
specimens  shown  in  Figure  54(c  and  d).  To  ensure  bond 
integrity,  the  metallic  shims  were  appropriatelv 
cleaned  and  coated  with  a  thin  laver  of  4617  epoxv 
before  insertion  into  the  laminate.  In  each 
specimen  two  16-mil  and  two  20-*;' 1  thick  shims 
were  used. 

Stiffened  Compression  Panel 


i  /..i.t  Ti?e  cutout  reinforcement  consisted  of  3  parts:  one  8- 
fiL  *±/5'  °'  ?'  9'  °'T45  )  which  was  bonded  to  the  skin  and  two 
stringer  cap  load  transfer  plates  of  a  4-plv  (+45?  *  45  ) 

configuration  .  See  Figure  55.  Seven  hat-shaped  stringers  were 
used  with  a  stringer  spacing  of  2.25  in.  The  legs  of  the  3 

ac2o^od^£in?hrS  ufre  "5°f9led"  in  design  and  fabrication  to 
accommodate  the  skin  reinforcement  plate.  Stringers  and 

reinforcements  were  bonded  using  Epon  921  adhesive  witt  a  room 

eu?nn£atUre  CUre  t0  miniinize  thermal  stresses.  The  shape  of  the 

structure?8  representative  of  “n  access  opening  in  an  airframe 


C.  Test  Results 


Photographs  of  the 
‘52(a),  56  and  57  with 
XIX. 


panels  after  test  are  shown  in 
a  summary  of  failure  loads  given 


Fiqures 
in  Table 
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Figure  55.  Stiffened  Cor  ossion  Panel  with 
Reinforced  Square  Cutout. 
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Figure  56.  Cutout  Panels  after  Test 


Figure  57.  Reinforced  Cutout  Panels  after  Test. 


TABLE  XIX 


COMPARISON  OF  FAILURE  LOADS  (LB/IN. ) 


Panel 

Predicted 

Measured 

CONTROL 

3240 

2420* 

UNREINFORCED 

1350 

1300 

BONDED  EXTERNAL  STEEL 

2200 

2250 

INTEGRAL  COMPOSITE  (1) 

2480 

2260 

INTEGRAL  COMPOSITE  (2) 

2480 

2200 

INTEGRAL  TITANIUM 

2920 

2020 

INTEGRAL  STAINLESS  STEEL 

2760 

1990 

STIFFENED  COMPRESSION 

1430** 

1610 

3900*** 

*  Failed  near  end  attachment. 

**  Southwell  prediction  for  critical  buckling  load. 
***  Compression  strength  prediction. 
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With  reference  to  control  panel  failure  at  75ft  of  predictor! 
strength  as  shown  in  Table  XIX,  it  must  be  noted  that  failure 
did  not  occur  in  the  gage  section;  sse  Figure  52(a).  Analysis 
had  shown  that  a  complex  biaxial  stress  state  exists  in  the 
corner  of  the  panel  near  the  end  attachments  as  a  result  of 
clamped  boundary  conditions.  The  analysis  in  this  region 
indicated  that  failure  should  initiate  at  83ft  of  ultimate 
laminate  tensile  strength.  This  analytical  value  is  in 
reasonable  agreement  with  the  experimental  value  of  75ft 
considering,  as  well,  that  nominal  material  properties  were  used 
previously. 

Considerable  strain  gage  instrumentation  was  used  on  the 
control  panel  as  seen  in  Figure  52(a).  All  load-strain  curves 
were  linear  to  failure  and  showed  excellent  strain  uniformity 
across  the  panel  midsection  with  negligible  bending.  Predicted 
and  measured  results  of  the  control  panel  test  are  summarized  in 
Table  XX. 


TABLE  XX 


CONTROL  PANEL  TEST  RESULTS 
"THORNEL"  50S/4617,  (+15,  90,  *15  ) 


Property 

Predicted 

Measured 

X 

(ksi) 

95.5 

69.5* 

T 

(  € 

)  (  in. /in.) 

4100 

2970 

1 

ult 

6 

E 

(10  psi) 

23.2 

23.4 

1 

V,  _ 

0.259 

0.222 

12 

V 

(ft) 

Same  - - 

-  57.4 

F 

V 

(ft) 

- 

0 

V 

♦Failure  near  end  attachment 
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,(fTI,  Ta^Je  XXI  summarizes  predicted  undirectional  data  for  the 
Thornel  50S/4617  system  and  complete  properties  for  the  test 
iaminate.  The  unidrectional  properties  have  been  corrected  as 
noted  based  on  tests  of  a  5-ply  unidirectional  panel.  The  latter 
data  are  presented  in  Table  XXJI  . 

Tne  unreinforced  cutout  panel  failed  at  54 %  of  the  load  that 
the  control  panel  carried  as  may  be  seen  ftom  Table  XIX. 

WaL  catast,rophic  and  sudden,  it  was  believed 
that  tensile  failure  first  took  place  at  the  left  side  of  the 

hoie  foliowed  by  crack  propagation  out  toward  the  free  left  edge. 

ffcn TifUre-56  ai*  Eccentrif  loading  and  bending  would  then  have 

resulted  m  failure  at  the  right  side  of  the  hole  which 
propagated  upward  and  to  the  right.  Amost  all  strain  gages  were 
increasing  load  up  to  1000  lb/in.  Above  this  load, 
nonlinearities  were  pronounced,  especially  in  the  vicinity  of  the 
hole.  Typical  strain  responses  are  shown  in  Figure  58.  .  Gage 

9  shown  therein  measured  the  highest  strain,  3700 fi  in. /in. 

It  can  also  be  seen  from  Table  XIX  that  the  reinforced 

hhe°?t  Hpanel?  .failed  at  no  less  than‘  82*  and  as  high  as  93%  of 
the  load  carried  by  the  control  panel.  Of  these  the  highest 

ailure  load  was  achieved  with  the  composite  (1)  design  As 

leZfci"  Higure  V  failure”  ocourred  fche  ends  of  the  fell?* 
reinforced  specimens  and  were  ~  30%  less  than  predicted  values. 

These  failure  modes  were  attributed  to  the  fact  that  the  metal 
shims,  although  0.016  in.  thick,  ended  abruptly  and  caused  local 
stress  concentrations  which  initiated  failure.  With  the  graphite 
reinforced  panels,  this  problem  was  minimized  by  curing  the 

transiti:onrtSp^f  panela  simultaneously,  thus  obtaining  smoother 

ifee  feictfefefe  graPhlte  reinforced  burred 

A  comparative  failure  load  and  reinforcement  weight  analysis 

to  thlTof  fLFigT  i9'  failUre  loads  have  bee^  normalized 

hat  of  the  control  panel  and  the  normalized  weights  were 

computed  as  total  reinforcement  weight  divided  by height  of 

co^Mal  ^m°Ved  the  CUt°ut*  This  comparison  shows  Sat  a 
d  b}®  Weig^  swings  can  be  achieved  while  maintaining 

o£  graphite  i„g 

--  « “y  sirs  %;isatj™Ltx°e"1 

delaminated  as  shown  in  Figure  60  before  panel  failure. P  The 
delammaiion  shown  extended  to  a  depth  of  ~0.5  in.  into  eSch  of 

feonifefe  T  7eI7  durlng  test  at  a  .S2.  Of 

43,000  psi  or  62%  of  ultimate.  Visual  inspection  revealed  a 
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TABLE  XXI 

PREDICTED  "THORNE  L"  50S/4617  PROPERTIES 


Property 

Unidirectional 

Laminate 

Test  Laminate 

5  (±15,  90,  =(=15°) 

VF 

(%) 

Same  -«i — 

- 57 

t 

(in.) 

- 

0.034 

Ej  ^ 

33.4 

23.2 

e2 

>  (10®  psi) 

1.0 

7.5 

Gl2  _ 

0.5 

2.1 

*12 

0.304 

0.259 

*21 

0.012 

0.084 

xT 

128  (2) 

95.5 

*c 

100 

75.2 

yt 

>  (ksi) 

2.5 

25.6 

YC 

28.0 

22.5 

T  J 

6.3 

21.2 

(1)  Ei 

(2)  XT 


1.13  x  rule  of  mixtures  value 

0.895  x  rule  of  mixtures  value  J  See  Tab*e  XXII 


TABLE  XXII 

PREDICTED  AND  MEASURED  “THORNE L”  50S/4617 
UNIDIRECTIONAL  PROPERTIES 


Property 

Predicted^) 

Measured 

XT 

(ksi) 

161.0 

144.0 

E1 

(10®  psi) 

32.0 

36.1 

vF 

<%) 

Same  — 

—  65 

VV 

<%) 

- 

0.5 

(1)  Based  on  rule  of  matures  with  Vp  »  65%,  cr  Fiber  =  245  ksi 
(yarn  strength)  and  Ep  -  48  x  106  psi  per  Union  Carbide  data 
this  lot. 


131  - 


132 


igure  58.  Typical  Strain  Responses  from  Cutout  Panel  Test 
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Thornel"  50S/4617  Specimens 
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Figure  6G„  Typical  Edge  Delamination 


failure  approximately  in  the  center  of  the  90°  ply  and  normal  to 
the  plates.  Strain  gage  data  did  not  show  evidence  of  this 
delamination.  Due  to  the  antisymmetry  of  the  +15°  plies,  the 
panel  edges  were  observed  to  curl  after  test  as  shown  in  Figure 
60. 


Ends  of  the  compression  panel  were  potted  to  simulate 
clamped  end  conditions  during  test.  The  failure  load  was  1610 
lb/in.  as  compared  to  3900  lb/in.  predicted  for  a  pure 
compressive  failure.  A  Southwell  plot  of  the  panel's  midpoint 
transverse  deflection  indicated  the  panel  to  he  buckling  critical 
at  a  load  of  1430  lb/in.  This  load  was  reasonably  close  to  the 
measured  failure  load.  Failure  did  rat  occur  in  the  reinforced 
region, .  rather  it  was  a  combination  shear/compressive  tvne 
traversing  across  the  panel  approximately  1  in.  above  the 
reinforced  region. 

D.  Test-Theory  Comparisons 

A  discrete  element  analysis  was  performed  for  each  of  the 
unstiffened  panels  tested.  Due  co  symmetry  of  the  plate  and 
loading,  only  one  quadrant  of  the  plate  was  analyzed,  and 
symmetric  boundary  conditions  were  imposed  along  two  axes. 

The  plate  was  idealized  as  a  two-dimensional  (plane  stress) 
assemblage  with  39  quadrilateral  plate  elements.  To  simulate 
uniaxial  tensile  test  conditions,  a  uniform  displacement  (ua)  was 
imposed  in  the  pulled  direction  while  restraining  displacements 
in  the  transverse  direction. 

A  margin  of  safety  analysis  was  performed  for  each  panel, 
and  the  predicted  ultimate  loads  were  calculated.  A  comparison 
of  the  predicted  and  measured  failure  loads  was  given  in  Table 
XIX  showing  fair-to-good  agreements.  Figure  6 l  shows  a 
comparison  of  predicted  and  measured  strain  response  for  the 
unreinforced  cutout  panel  showing  excellent  aaree’Tnt  between 
discrete  element  and  measured  strains.  Similar  agreements  were 
obtained  for  the  other  panels.  Analytically  predicted  stress  and 
strain  concentrations  at  the  edge  of  an  unreinto-ced  circular 
cutout  on  an  infinite  plate  are  shown  in  Fioure  62..  These 
concentrations  ^are  equal  with  a  value  of  4.75  at  9  =  90°  and 

T*no\at  for  the  test  Panel  orthotrop'/  ratio  (Ei  /E?  = 
3.J9).  Elsewhere  around  the  hole,  stress  and  strain  "K"  factors 
are  unequal  in  contrast  to  those  for  the  isotropic  material. 
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